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ABSTRACT 



The problem of finding the shortest path between two 
points in a plane containing obstacles is considered. The set 
of such paths is uncountably infinite, making an exhaustive 
search impossible. This difficulty is overcome by reducing 
the size of the search space. The search is first restricted 
to a countably infinite set by focusing attention on the set 
of homotopy classes. By applying simple optimality 
principles, a finite list of such classes is obtained whose 
union contains the shortest path. This process of 
simplification is accomplished by modeling the topology of the 
region with a graph. Optimality principles come into play 
during a graph traversal which is used to produce the finite 
search list. In addition, a computational investigation of 
two methods by which homotopy classes can be named is 
discussed, and properties of the graph models are 
investigated. The thesis of CPT Andre M. Cuerington, U.S. 
Army, calculates the actual shortest path using the search 
list produced here. 
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THESIS DISCLAIMER 

The reader is cautioned that computer programs developed 
in this research may not have been exercised for all cases of 
interest. While every effort has been made, within the time 
available, to ensure that the programs are free of 
computational and logic errors, they cannot be considered 
validated. Any application of these programs without 
additional verification is at the risk of the user. 
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INTRODUCTION 



I . 

The industrial conununity relies today, in increasing 
measure, upon robots to handle a multitude of tasks. One goal 
is to develop robots, and fixed robot arms, with the ability 
to roam freely among obstacles. In this regard, the 
calculation of shortest paths is of obvious importance. 

A . THE PROBLEM 

As robot arms move, they must vary their configurations to 
position their end-effector or "hand". Once obstacles are 
introduced, the freedom of movement of the manipulator may 
become drastically reduced. The problem considered here is 
one of finding the shortest path from a starting point to a 
destination point in a planar region containing obstacles. 
This paper addresses several parts of the solution to this 
problem. 

B. OVERVIEW OF THE SOLUTION 

Between any two points in the plane, a path joining them 
may be chosen from an uncountably infinite set of 
alternatives. It is desirable to choose the shortest path 
from this set which, due to the presence of obstacles, may not 
be a straight line. An exhaustive search of the collection of 
possible paths is impossible, so another method must be 
developed. 
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The set of possible paths is partitioned into a collection 
of equivalence cl asses--mutual 1 y exclusive subsets which 
collectively exhaust the partitioned set. This partition is 
induced when an equivalence relation is defined on the set of 
possible paths. This set of equivalence classes is countably 
infinite . 

Naming conventions are established which associate with 
each class a character string which allows the classes to be 
represented in a computer. Next, a finite list of candidate 
classes is produced to search. 

A heuristic is then applied to this finite list of classes 
to place them in an order which saves computational effort. 

The final step begins with this ordered list of candidate 
classes and solves the shortest path problem, class by class, 
in listed order. Savings of computational effort are realized 
when the class-by-class solution process can terminate without 
exhausting the ordered list of candidate classes. 

The key ideas in the above overview are now considered in 
somewhat greater detail. 

C . THE APPROACH 

1. The Topology 

Let P be the uncountably infinite set of all 
continuous, obstacle-avoiding paths from a starting point a to 
a destination point z. If Pj and pj are coterminal paths in 
P, that is, paths which have the same starting point and 
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destination point, we say that is homotopic to Pj if p^ can 
be mapped to Pj under a continuous function (with both 
endpoints fixed in place) without encroaching on any obstacles 
[Ref. l:p. 223]. The homotopy relation is reflexive, 
symmetric and transitive and therefore defines an equivalence 
relation [Ref. l:p. 223]. This relation induces a partition 
of P into a countably infinite collection of equivalence 
classes, known as homotopy classes [Ref. l:p. 223]. 

2. Establishing the Naming Convention 

In order to name homotopy classes, a reference frame 
is established to represent the topological relationships in 
the plane with obstacles. For a given path, p, a string of 
characters, R(p), is recorded which encodes information 
concerning the relationship of the path to obstacles. 

Two algorithms are then presented which accept R(p) as 
input. These algorithms have the following important property 
which implicitly defines names for the homotopy classes: if 
P^ and Pj are coterminal paths and R(Pj) and R(pj) are input to 
either of the two algorithms, then the outputs are identical 
if and only if Pj is homotopic to Pj [Ref. 2]. 

Appendix A presents a computational investigation 
which closely examines the methods employed by these two 
algorithms in producing the names for these homotopy classes. 

3. Generation of Names for Candidate Equivalence Classes 

An edge-labeled graph is constructed which models 

topological relationships within the region, where the nodes 
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represent subregions induced by the reference frame. Nodes 
are connected if their corresponding subregions are adjacent. 
The names of the desired homotopy classes may be obtained by 
traversing this graph. This paper presents algorithms both to 
construct the graph and to traverse it. This traversal 
produces a list of candidate classes, which contain all paths 
of minimal length. 

The location of the origin affects the manner in which 
the plane is divided into wedges. In this regard, the graph 
which is created, given the location of the origin in one 
instance, may not be the same as the graph obtained if the 
origin is moved to a new point. This thesis presents a method 
to count the number of possible different graphs, given n 
obstacles, and a method to construct them. 

4. Heuristic Ordering of Candidate Homotopy Classes 

For each class on the list, a lower bound on the 
length of its shortest representative path is constructed. 
The list of classes is then arranged into increasing order of 
these bounds. 

To obtain the bounds, a point is first chosen within 
each obstacle. A contraction deformation is then applied to 
each obstacle to "shrink” the obstacle so that the chosen 
point represents the obstacle. 

Then some class is fixed and its shortest path is 
examined. As all obstacles are simultaneously contracted to 
their representing points, this shortest path has a limiting 
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position which is polygonal. The length of this polygonal 
path is the lower bound associated with that class. The 
length can be calculated from the class name and representing 
points without explicitly defining any contraction 
deformation . 

5. Class by Class Solution of the Shortest Path Problem 
In the final step of the solution, the classes on the 
ordered list of candidates are considered. The first class is 
removed from the list, i.e., that class with the smallest 
associated bound. The true shortest path in the named class 
is found by reversing the contraction previously applied to 
the obstacles, thereby transforming the polygonal path whose 
length provided the lower bound into the true shortest path. 
If the length of this path is smaller than the bound 
associated with the class on the top of the remaining list, 
the search is stopped. Otherwise, the first class is removed 
from that remaining list and the procedure is repeated. This 
procedure continues until the condition specified above is 
met, and the shortest path has been found. 

D . SUMMARY 

The solution to the problem of searching for the shortest 
path between two points in the plane with obstacles begins 
with consideration of a set of paths which is uncountably 
infinite. Through the homotopy relation, this set is reduced 
to a countably infinite set. The new set is further reduced 
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to a finite list by modeling the region with a graph and 
applying an optimality principle. This final list contains 
only those homotopy classes containing paths which are not 
sel f -crossing . From this list a shortest path is found. 

Computational effort is further reduced through the use of 
a heuristic which orders the list of candidate classes by 
increasing order of their lower bounds. The heuristic used 
also facilitates a methodical search of the classes while 
solving the shortest path problem. The use of this heuristic 
does not, however, imply that the solution is approximate. 
The solution to the problem will be exact using this method. 

E. THE CONTRIBUTION 

Chapter II of this thesis presents the methods used to 
establish a reference frame given the plane with obstacles and 
an algorithm which is used to generate the homotopy class name 
of a given path. 

Chapter III introduces a graph which is used to model the 
region. An algorithm is then presented which is used to 
traverse this graph to produce the list of candidate homotopy 
classes. Using the results of this algorithm, the thesis 
submitted by CPT Andre M. Cuerington, U.S. Army [Ref. 3], 
shows how to find the solution to the shortest path problem. 

Chapter IV introduces some properties of the graph which 
are used to model the region. 
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chapter V employs counting methods to determine the number 
of non-isomorphic graphs which may be encountered given a 
plane with some number of fixed obstacles. The number of non- 
isomorphic graphs represents the number of truly different 
reference frames which may be constructed in the region, thus 
producing different sets of names for the homotopy classes. 

Finally, Appendix A presents a computational investigation 
which addresses the question of whether or not two algorithms 
used to name homotopy classes actually produce the same 
resul ts . 
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II. NAMING EQUIVALENCE CLASSES 



A. INTRODUCTION 

In this chapter, a labeling scheme is established to 
represent the equivalence classes. This notation is used 
throughout this paper to organize the computational search for 
the shortest path. 

An algorithm is presented that from a character 
representation of a path, p, names the path's respective 
equivalence class. It can be shown that, after being 
processed by the algorithm, different character 
representations of paths Pj and Pj yield the same output 
exactly when the two paths are in the same equivalence 
class [Ref. 2]. Thus, the name of a class will be the string 
obtained by applying the algorithm to any representative of 
the class. 

The computation of class names depends on a reference 
frame, which in turn depends on a collection of obstacles. 
Although more than one reference frame can be drawn for a 
particular collection of obstacles, the choice of a particular 
frame fixes the representation of all homotopy classes. 

Once a reference frame is developed, the name of the 
homotopy class for a path can be determined by a two-step 
procedure. First, a character string R(p) is calculated which 
encodes certain information about the path taken through the 
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obstacles. Second, this character string is accepted as input 
to an algorithm, which then produces a name for the 
equivalence class to which the particular path belongs. This 
algorithm produces the class names in terms of the same 
alphabet used to create the initial character string. The 
results from the algorithm described. Algorithm 1, will be 
used throughout this analysis as the class name associated 
with a given path. 

B. CONSTRUCTION OF REFERENCE FRAME 

Let bjj be an arbitrary point in obstacle Bj^, k = 1,2,. ..,n, 
where n is the number of obstacles in the region. A point c 
is chosen and a line drawn through each b|^, infinite in extent 
in each direction, and having the following properties: there 
is an open disk, 6, centered at c such that SHBj^ = 0 for all 
k = l,...,n, and the n lines connecting c with the points b]^ 
are distinct. Such a point c can always be found as the above 
two conditions are satisfied by any point in the planar region 
that is neither on an obstacle nor on the n(n-l)/2 lines 
determined by pairwise choices of distinct b|j. 

To draw a reference frame, n lines are first constructed 
joining c with each b|^. The line from c to bj^ is labeled as 
L|j. Each line is then partitioned into two directed, semi- 
infinite rays and one finite length line segment. The ray 
emanating from c in the direction away from b|^ is called . 
The ray emanating from bij and away from c is denoted Cjj . The 
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remaining line segment, [c,bjj], is also denoted Qj^. The 
reference frame is the collection of line segments and rays so 
constructed, as illustrated in Figure 2.1. 




Figure 2.1 A Reference Frame With n = 2 Obstacles 
C. RAW CHARACTER STRINGS 

A reference frame is fixed in a region T. Let a and z be 
points in T, and p be a directed path in T from a to z. Then 
the raw string of p, denoted R(p) , is defined to be the 
ordered sequence of characters obtained by following p from a 
to z and recording, in order, the names of the rays that are 
crossed . 

Two special cases must be addressed to make the above 
definition complete. First, in the case that p crosses no 
reference rays, we set R(p) = e, where e denotes the empty 
string. Second, if p crosses through c (simultaneously 
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crossing all Qj , the names of all Qj will be recorded in order 
of increasing subscript. 

Thus, raw strings are of the form 

R(p) = XjXj. . .x„ 

where each Xj belongs to the alphabet 

A = {e, Qj, Qj , . . . , ct„. Bp..., flj} . 

The character x^ (j = 1, 2,..., m) is the name of the 

reference ray crossed by p. Figure 2.2 shows a pair of paths 
connecting the two points a and z in the region with two 
obstacl es . 




Figure 2.2 Paths Q 2 ^i 

D. ALGORITHM 1 

Algorithm 1 accepts as input a raw string R(p) for any 
path p in T and any reference frame. The output is a string, 
denoted C(R(p)), of characters also chosen from the 



11 



alphabet A. This output string C(R(p)) is the name of the 
homotopy class to which p belongs. 

The algorithm is presented in terms of two functions. The 
first is the sorting function o. Let S = Xj . . . be any string 
over A. If S contains a two character substring XjjXjj^j = 
with i<j, then o(S) is the string which results by reversing 
the order of the leftmost such substring. Figure 2.3 depicts 
two such strings and makes clear that such paths are 
homotopic. If S contains no such two character substring, 
then o ( S ) = S . 




Figure 2.3 Homotopic Alpha Strings Q 2 Qj and 0^02 

Repeated application of o sorts all substrings that 
consist entirely of into non-decreasing order of 
subscript. Z(S) is defined to be the string which results 
when this sorting is complete. 
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The second function is the cancellation function %• If 
string S contains a two character substring with X|j = Xj^^j 
then x(S) is the string which results by removing the leftmost 
such two character substring. Otherwise x(S) = S. Such 
cancellation is intuitive, for if a reference ray is crossed 
twice consecutively, then that is equivalent to not crossing 
at all. Thus repeated application of % cancels all pairs of 
adjacent like occurences of identical characters. Let X(S) be 
defined as the string which results when all such possible 
cancellation is complete. 

With these definitions complete. Algorithm 1 is given in 
Figure 2.4. 



BEGIN 
k - 0 
S|t - R(p) 

WHILE S|^ not equal to X(i:(S|^)) 

k k+1 

S)j - X(Z(S]j.l)) 

END WHILE 
C(R(p)) - S), 

END 



Figure 2.4 Algorithm 1 



The output from Algorithm 1, C(R(p)), is called the 
canonical representation of p and is the unique name for the 
homotopy class of p [Ref. 2]. That name is used throughout 
this paper. 
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Ill . 



A GRAPH MODEL FOR NAMING CLASSES 



This chapter presents a graph used to model the 
topological relationships in the region. The graph is then 
used to produce a finite list of names of candidate homotopy 
classes for subsequent searching. 

r 

A. MODELING THE TOPOLOGICAL SPACE 

The reference frame, once established, divides the region 
into 2n "wedges". Each wedge is represented by a node. Two 
nodes are connected by one edge if the wedges they represent 
share a boundary of positive length corresponding to only an 
a ray, or two edges if the two wedges share a boundary 
corresponding to both a and B rays. No edge connects two 
nodes if the wedges they represent have no common boundary 
other than the origin, c. Figure 3.1 illustrates this 
procedure of replacing the wedges by their associated graph 
with n = 2. 

Figure 3.1 Establishing the Nodes and Edges of the Graph 
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Each edge of the graph is then labeled with an or a A); • 
An Ojj is assigned to an edge if the incident nodes represent 
wedges sharing a boundary created by an 0 |j reference line. 
The labels are similarly assigned. When parallel edges 
connect adjacent nodes the outermost edge is labeled with , 
the innermost edge with ajj . This labeling coincides with the 
fact that the ray for each obstacle is on the outboard side 
of the obstacle with respect to the origin, c. Figure 3.2 
illustrates the edge labeling for n = 2. 





Figure 3.2 Edge Labeling for Graph with n = 2 

In constructing the model , wedges and nodes are numbered 
counterclockwise using upper case Roman numerals. 

We define a ring graph to be such a graph used to model a 
region. Figure 3.3 illustrates the region with n = 2 and its 
associated ring graph. 
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If a and z are points in the topological space, T, and p 
is any path from a to z, then the string of characters, R(p), 
which represents that path can be produced by tracing p from 
a to z in T and recording R(p) as the list of rays crossed, as 
defined in Chapter II. However, R(p) can also be produced by 
traversing the ring graph and visiting the nodes which 
correspond to the wedges visited by p in T and recording the 
name of each edge as it is traversed. Figure 3.4 shows a path 
in a region where n = 2, and illustrates the path in the 
corresponding ring graph. 

In constructing the model , it becomes unnecessary to 
generate a complete graph on the 2n nodes by eliminating the 
origin as a boundary between any two wedges. This may be 
justified by observing that any path passing directly through 
the origin could be perturbed in such a manner that any 
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analysis of the path would remain unchanged. This situation 
is illustrated in Figure 3.5. 




R(p) - 0^2 P 1 

Figure 3.4 Path Represented by Ring Graph for n = 2 




Figure 3.5 Perturbation of a Path Through the Origin 
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B. PRODUCING CLASS NAMES USING GRAPH TRAVERSAL 



Any path p in T has a character string representation 
associated with it which may be obtained by traversing the 
graph. A list of character strings which name the homotopy 
classes is desired. This list of names is finite and includes 
only the names of those classes whose representative of 
shortest length is not sel f -crossing and starts and finishes 
in the appropriate wedges. 

A list of character strings is constructed by beginning at 
some node and searching the graph using a modified Breadth 
First Search (BFS) algorithm. This search utilizes a first- 
in/first-out queue which stores the current list of character 
strings and an element which denotes the node in which each 
string terminates. 

To illustrate the procedure, we begin with an empty string 
and the starting node of the graph on the queue. Then we copy 
the first string and the node denoting the current position 
from the queue. A new character which represents an edge 
which may be traversed in exiting the current node is added to 
the string. This new string is both printed and added to the 
bottom of the queue. All similar strings are generated by 
traversing each of the edges which may be traversed from the 
current node. The procedure is then repeated for the next 
string on the queue. Figure 3.6 presents this search in 
algorithmic language. 
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k := bottom of queue 
Q := top of queue 
v() := vertex on queue 
s() := string on queue 

BEGIN 
k - 1 
Q - 1 

v(k) •- start vertex 
s(k) - e 
WHILE (Q<k) DO 

string *- s(Q) 
vertex ♦- v(Q) 

FOR EACH EDGE OF THE FORM (vertex, vertex2) DO 
k •- k + 1 

char «- label [( vertex , vertex2)] 
s(k) *- string:char 
v(k) •- vertex2 
OUTPUT s(k) , v(k) 

Q - Q + 1 
END WHILE 
END 

Figure 3.6 Modified BFS Algorithm 



In the absence of some criterion for terminating this 
search, an ever-growing list would be produced. A new string 
would be produced for every node which is adjacent to the 
current node. It is easily seen that this infinite list would 
eventually include the character string representation of 
every possible path in T which starts in some given wedge. In 
order to accomplish the task of reducing this list to a finite 
list of names of candidate homotopy classes, some stopping 
criterion is required. 

By examining the list of strings, it can be seen that some 
of the strings are not desirable because they are 
representative of classes which contain only paths which are 
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sel f -crossing or are subject to cancellation under 
Algorithm 1. These paths should be removed from the list. 
Since all of the extensions of these excluded strings are also 
undesirable, they too should be excluded. The algorithm then 
terminates after producing a finite list. That the list is 
finite is substantiated by using the Pumping Lemma of finite 
state machines. This lemma states that for any language, an 
infinite string is possible if and only if cycles are 
permitted [Ref. 4:p. 77]. The language here consists of the 
R(p) strings. A cycle in the language is representative of a 
loop around any obstacle or set of obstacles. Therefore, if 
we remove the cycles, a language with a finite number of words 
is produced, which in this case is a finite number of paths. 

Two such undesirable strings are shown in Figure 3.7. 
These strings represent • paths which either wrap around an 
obstacle or cross themselves by circling two obstacles in 
figure eight fashion. 
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By performing several checks on each string as it is 
generated by the search algorithm, those strings which have an 
undesirable form may be "sieved" from the list. In other 
words, these strings will not be placed on the queue, and will 
therefore no longer be considered in the search. The criteria 
which are applied to each string to accomplish this sieving 
are given below. 

The first criterion is that any string containing a 
substring of the form a]j,Gjj,a]j or 13]^ , O), , 13]^ is not placed on the 
list. Such a string represents a class whose shortest 
representative wraps around an obstacle and is self crossing. 
The second criterion is that strings containing substrings of 
the form or Bj^ , , Bjj , represent classes whose 
shortest representative is also self crossing. They too are 
not placed on the list. Figure 3.8 provides examples of paths 
which produce substrings of these types, and also demonstrate 
how the shortest representative in each class crosses itself. 

A third criterion eliminates from the list those strings 
which contain two like a's not separated by any B. Such 
strings are subject to cancellation as defined in Algorithm 1 
presented in Chapter II. Once cancellation takes place the 
resulting string merely duplicates some string already on the 
list. 

Another criterion introduces the restriction that the 
algorithm may not traverse an edge of the graph two times in 
immediate succession. Therefore, generation and cancellation 
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of such strings is not required. Also, duplication of strings 
can not occur. 

Finally, those strings which do not terminate at the 
desired node are eliminated from the list. 

What remains on the list upon termination of the search 
given these sieving rules are the names of only those 
candidate homotopy classes which are to be used in the 
subsequent search for £ shortest path. 



a 





Figure 3.8 Examples of Self Crossing Representatives 
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IV. PROPERTIES OF THE RING GRAPH 



This chapter introduces some of the properties of the ring 
graphs obtained when the number of obstacles is greater than 
one. The case in which there is a single obstacle is 
addressed briefly at the end of the chapter. 

As stated in Chapter II, each reference line, Lj^, in the 
region containing n fixed obstacles is partitioned into three 
parts--two semi-infinite rays and one finite line segment. 
Since there is a one-to-one correspondence between these 
components and the edges of the ring graph, it is easily seen 
that the total number of edges in the graph will equal 3n for 
any n > 2 . 

Now, we superimpose a ring graph over its corresponding 
reference frame, as illustrated in Figure 4.1. We first 
observe that, by definition, each Lj^ contains exactly one 
point bjj. That point contributes to the existence of two 
edges connecting the adjacent wedges. Likewise, since no such 
b]j lies on Lj^ on the opposite side of c, only one edge connects 
the adjacent wedges. Therefore, it is easily seen that in 
every ring graph with n > 2, there are always two parallel 
edges opposite a single edge. 
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Figure 4.1 Region and Ring Graph Superimposed for n = 3 

The number of edges incident with a node is called the 
degree of the node [Ref. 5:p. 86]. Consider the degree of the 
nodes of the ring graph. Every node in a given ring graph is 
of degree two, three or four. The sum of the degrees of 
the nodes in a graph is equal to twice the number of edges in 
that graph [Ref. 5:p. 86]. Since it has already been stated 
that the number of edges in the ring graph equals 3n, then for 
each ring graph, 

2n 

^ degreeivertex i) = 6n . (4.1) 
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since we are interested in the degree of incidence of 
individual nodes we rewrite (4.1) in the form 

4 

6n = ^ iinumbez of vertices of degree i) . (4.2) 

j=2 

It is an easy combinatorial result that in any graph, 
there are an even number of nodes of odd degree 
[Ref. 5:p. 91]. Therefore, in any ring graph, the number of 
nodes of degree three is even. 

Define g, h and j to be the number of nodes of degree two, 
three and four, respectively. Then the total number of nodes 
is 

g + h + j = 2n , (4.3) 

and (4.2) may be rewritten as 

2g + 2h + 4j = 6n . (4.4) 



Multiplying both sides of (4.3) by 3 and combining equations 
(4.3) and (4.4) yields 

3 ig + h + j) = 2g + 3h + 4j , (4.5) 

The solution to (4.5) shows that g = j or simply that the 
number of nodes of degree two equals the number of nodes of 
degree four in any ring graph with n > 2. It may also be 
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verified by examining any ring graph that every node of degree 
two is directly opposite a node of degree four. Hence, there 
must be an equal number of each. 

All cycles in ring graphs are on an even number of nodes. 
Therefore these graphs may be classified as being of bipartite 
type. That is, the set of nodes may be partitioned into two 
sets such that all edges of the graph connect a node in one 
set to a node in the other set [Ref 5:p. 240]. Since the 
nodes are labeled sequentially in a counterclockwise 
direction, the two sets in this partition consist of those 
nodes whose labels are even integers and those whose labels 
are odd integers, respectively. 

This bipartite structure allows certain observations 
concerning the distance travelled in moving from node to node 
through the ring graph. ' Assigning an equal weight of one to 
each edge, the distance travelled along a path is simply the 
total number of edges traversed (with the multiplicity of 
traversals of individual edges included). It is easily seen 
that any path which begins and ends at nodes whose labels are 
of like parity will have even length. On the other hand, a 
path which begins and ends at nodes whose labels are of 
differing parity will have odd length. 

Consider a chain of m nodes of degree three. For m odd, 
the chain appears as that shown in Figure 4.2(a). For m even, 
the chain takes on one of the forms shown in Figure 4.2(b). 
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0^-0 



< 03 — <0 



o — O — <0 

(a) (b) 

Figure 4.2 Chains of Nodes of Degree Three 

Consider a chain such as the one illustrated in 
Figure 4.2(a). It is easily seen that in concatenating this 
chain with nodes of degree two and four, it is possible to do 
so with only a degree two node on one end and only a degree 
four node on the other. Figure 4.3 illustrates the extension 
of such a chain. 






Figure 4.3 Concatenation of Even Degree Nodes 

Using a similar construction, chains such as those 
illustrated in Figure 4.2(b) may be used to show, as seen in 
Figure 4.4, that only nodes of the same degree (two or four) 
may be concatenated on both ends. 
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^ — cOd — e 



oo — o — cc>0 

Figure 4.4 Concatenation of Even Degree Nodes 

These facts lead to the observation that any string of 
nodes of degree two and any string of nodes of degree four 
must be separated by an odd number of nodes of degree three. 
Furthermore, if two nodes of degree two (or four) are 
separated by nodes of degree three, the number of degree three 
nodes must be even. 

It is possible to represent a given ring graph by an 
adjacency matrix, A. The resulting matrix is a symmetric 
2n X 2n matrix in which the (i,j) element is a 1 if there is 
a single edge connecting nodes vj and Vj , a 2 if two parallel 
edges connect the nodes, or a 0 if the nodes are not adjacent. 
Figure 4.5 shows a ring graph with n = 3 and its associated 
adjacency matrix. By virtue of the conventions used to label 
the nodes, it is observed that all nonzero entries in the 
adjacency matrix are on the superdiagonal , the subdiagonal and 
in the (2n,l) and (l,2n) positions. (All entries on the main 
diagonal of A are zero since no ring graph contains any edges 
which connect a node to itself.) Hence, the matrix becomes 
very sparse as n becomes large. 
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A basic property of an adjacency matrix is that the sum of 
the values in the i^^ row (column) of the matrix is equal to 
the degree of node Vj in the graph [Ref. 5:p. 88]. 

As stated above, the conventions used to label the nodes 
of the ring graph indicate that any two nodes whose labels are 
of like parity must be an even distance apart. By the 
bipartite property of the graph, all (i, j) entries in A, where 
i and j are of like parity will necessarily be zero. 




0 0 1 
0 0 0 
2 0 0 
0 1 0 
1 0 2 
0 2 0 



Figure 4.5 A Ring Graph and its Adjacency Matrix, A 



Another theorem from graph theory states that the number 
of paths of length m from node vj to node Vj is the (i,j) entry 
of the m^^ power of the adjacency matrix. A*" [Ref. 5:p. 144]. 
This theorem leads to the further conclusion that when m is 
odd, all (i,j) entries where i and j are of like parity will 
be zero. On the other hand, when m is even, all (i, j) entries 
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for i and j of opposite parity will be zero. Naturally, all 
entries not mentioned above may or may not be zero, and depend 
on the actual layout of the region and its corresponding ring 
graph . 

In the case where there is only one obstacle, n = 1, a 
special case of the graph arises. Since one obstacle divides 
the region into two half planes, there are only two nodes in 
the graph with two parallel edges connecting them. Any two 
points in such a region may be connected by a path crossing 
either one a or one B edge. Therefore, this case is not 
considered in the analyses discussed in this paper. 
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V. 



NON -ISOMORPHIC RING GRAPHS 



A. INTRODUCTION 

The geometry of a region is fixed by the location of the 
obstacles. However, by varying the location of the origin, 
different reference frames may be established. These 
differences may lead to different ring graphs and therefore 
different homotopy class names. Figure 5.1 illustrates this 
possibility for n = 3, where two different reference frames 
result in different ring graphs and different class names for 
the same path. This leads to the following question: How 
many distinct ring graphs exist for a given collection of n 
fixed obstacles? 

It is shown in the sequel that the answer to the question 
may be found by determining the number of non-equivalent edge 
two-colorings of a 2n-gon subject to the restriction that 
opposite edges of the polygon are colored differently. 
Burnside's Lemma is used to determine this number. 

It must be pointed out that coloring a graph in the 
context of this paper is not the traditionally accepted 
definition of a graph coloring. Here, adjacent edges (or 
later, vertices) may be colored the same color. 

Once the number of different ring graphs has been 
determined, Polya's Pattern Inventory Theorem is used to 
construct them. 
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A complete example illustrating the procedures presented 
in this chapter is provided in Appendix C. Reference should 
be made to this example as new concepts are introduced. 







Figure 5.1 Effect of Reference Frame on Ring 
Graph and Class Names 

B. THE RING GRAPH/POLYGON CONNECTION 

"A graph is completely determined by specifying its node 
and edge sets. "[Ref . 6:p. 3] Any two graphs are isomorphic if 
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there exists a one-to-one correspondence between their node 
sets which preserves adjacency [Ref. 6 : p . 4]. The question at 
hand then becomes: How many non-isomorphic ring graphs exist 
for a given collection of n fixed obstacles? 

Each ring graph associated with a problem containing n 
obstacles may be modeled using a 2n-gon whose vertices 
correspond to the nodes of the graph. The edges of the 
polygon are then colored as follows: a black edge (bold line) 
connects those vertices whose corresponding nodes are 
connected by two parallel edges; a white edge (fine line) 
connects those vertices whose corresponding nodes are 
connected by a single edge; and no edge connects those 
vertices which correspond to non-ad jacent nodes. From this 
point forward, only single and double edges will be 
considered . 

It has been established that, in a ring graph, each set of 
parallel edges is opposite a single edge. Accordingly, the 
corresponding restriction on the polygon is that opposite 
edges must be colored differently. Figure 5.2 shows a ring 
graph for n = 2 and its associated 2n-gon. 

Since the configuration of a ring graph is defined by its 
node and edge sets, the location and/or spacing of its nodes 
is irrelevant. In constructing a graph then, the nodes may be 
thought of as being equally spaced around the ring. Therefore 
the 2n-gon associated with any ring graph may take the form of 
a regular polygon on 2n vertices. 
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Figure 5.2 Ring Graph and 2n-gon for n = 2 



With these ideas in mind, the problem of determining the 
number of non-isomorphic ring graphs given n fixed obstacles 
is the same as determining the number of non-equivalent edge 
two-colorings of a regular 2n-gon subject to the restriction 
that opposite edges are differently colored. 

C. COUNTING NON-EQUIVALENT EDGE TWO-COLORINGS 
1. Establishing the Permutations 

By coloring one edge of the 2n-gon, the restrictions 
which have been established previously automatically fix the 
color of the opposite edge. Therefore, of the 2n edges, n of 
them may be colored independently leading to a total of 2" 
possible colorings. 

Any pair of two-colorings of the polygon are 
considered equivalent if there is a rigid motion of the 2n-gon 
that maps the polygon onto itself [Ref 7:p. 335], and which 
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permutes one coloring into the other. Permissible rigid 
motions, or permutations, will be limited to rotations about 
the center of the polygon, reflections about diagonals drawn 
through opposite vertices, and reflections about bisectors 
drawn through opposite edges. Figure 5.3 illustrates each of 
these permutations using the square generated by n = 2. 
Examination of any 2n-gon reveals that there are a total of 2n 
rotations, n reflectiohs about diagonals and n reflections 
about bisectors. If G is defined to be the set of all 
permissible permutations on the set of edges of the 2n-gon, 
then it follows that |g| = 4n. 



\T 

\ ^ 

• — ^ ^ 




Figure 5.3 Permissible Rigid Motions of a Polygon 



2. Fixed Colorings Under Given Permutations 

Any permutation may be factored into a product of 
disjoint cycles [Ref 8:p. 69]. In order for a coloring of the 
polygon to be preserved under that permutation, all edges in 
a cycle must have the same color. Since opposite edges are to 
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be colored differently, they must be in disjoint cycles. Then 
all edges which are opposite the edges of a given cycle must 
themselves be in a disjoint cycle, which will be called the 
companion cycle. In order for a coloring to be fixed under a 
given permutation, each cycle must have a unique companion 
cycle, hence the number of disjoint cycles must be even. In 
such cases the edges in one cycle are colored black while 
those in its companion ' eye 1 e are colored white. Figure 5.4 
provides an example of a coloring which is fixed under a 
permutation, and shows the companion cycles for that 
permutation . 




Figure 5.4 Disjoint Cycles Produced by a Permutation 



We consider each of the three types of permutations. 
Throughout the following discussion 7(x) is defined to be the 
number of edge two-colorings which remain fixed under a given 
permutation, %. For simplicity, the edges of each 2n-gon will 
be numbered clockwise beginning at some fixed edge. 
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Rotations 



a . 

The identity permutation. Tig, is that motion which 
does not move the polygon. Hence, TCTig) = 2". The remaining 
rotations are through the first 2n-l multiples of (360/2n)° in 
the clockwise direction. 

Define a click to be (l/2n)^*’ of a full rotation 
of the polygon. A rotation, ti^, i = 1 ,...,( 2n- 1 ) , may be 
thought of as turning the polygon i clicks in the clockwise 
direction, thus moving each edge i positions along from its 
starting point. 

For the disjoint cycles generated by this 
permutation to have their coloring preserved, the number of 
disjoint cycles must be even. The edges in half of the cycles 
may be colored arbitrarily, while fixing the color of those in 
the companion cycles to the opposite color. Hence the number 
of choices to be made in coloring the edges of the polygon is 
one half the number of cycles, and f(Ti’) = If the 
number of disjoint cycles is odd then T(tc^) = 0. 

b. Reflections About Bisectors 

Once a bisector is drawn through two opposite 
edges of the polygon, the permutation for a reflection about 
that bisector yields a product of disjoint cycles consisting 
of two 1-cycles and (n-1) 2-cycles. Such a permutation is 
illustrated in Figure 5.5. In order for each cycle to have a 
companion cycle, (n-1) must be even. Hence, n must be odd. 
In this case coloring half the 2-cycles fixes the color of the 
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remaining companion cycles, and the number of possible choices 
for coloring the (n-1) cycles is (n-l)/2. One additional 
choice may be made in coloring those edges which lie on the 
bisector, resulting in a total of c = (n+l)/2 choices for 
coloring. Hence, when n is odd, 7(7t^) = 2'^, where 
i = 2n , . . . , ( 3n-l ) . 




Figure 5.5 Disjoint Cycles for Reflection About a Bisector 



For n even, the number of 2-cycles produced by the 
permutation is odd. This results in a cycle which fails to 
have a companion cycle. Hence no. coloring may be preserved 
under the permutation. Thus, when n is even, = 0, 
i = 2n , . . . , ( 3n-l ) . 

c. Reflections About Diagonals 

The permutation resulting from a reflection about 
an axis drawn through opposite vertices of the 2n-gon produces 
a product of disjoint cycles composed of n 2-cycles. 
Figure 5.6 illustrates such a reflection. It follows 
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immediately that when n is even, each cycle will have a 
companion and a total of n/ 2 choices may be made in coloring 
the 2n-gon. Then = 2°/^, i = 3n 4n-l ) . 

On the other hand, when n is odd, there will exist 
a cycle which does not have a companion cycle, and no coloring 
may be preserved under the permutation. This implies that 
when n is odd, = 0, i = 3n 4n-l ) . 




Figure 5.6 Disjoint Cycles for Reflection About a Diagonal 



3. Burnside's Lemma 

At this point, Burnside’s Lemma is used to calculate 
the total number, N, of non-equivalent edge two-colorings of 
a 2n-gon given the restriction that opposite edges must be 
colored differently. The lemma is as follows. 

Let E be the set of edges of a regular polygon and let 
G be a group of permutations, on the set E. Further, 
let C be any collection of colorings of E that is closed 
under G. Then the number, N, of non-equivalent colorings 
of the edges of the polygon is given by 
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(5.1) 



where |g| is the number of permutations in G and ¥(Jt) is 
the number of colorings in C which are left fixed by 
X . [Ref . 7 : p . 343 ] 

Given the information derived earlier, this lemma may 
be invoked to produce N. This number is given by 

4n-l 



Hence, the total number of non-isomorphic ring graphs given n 
fixed obstacles is obtained. 

D. PRODUCING NON- ISOMORPHIC RING GRAPHS 

Given the number of non-isomorphic ring graphs, N, 
obtained by applying Burnside's Lemma, the next step is to 
produce them. 

1. Using Partitions of the Integer n 



summands, mj,...,m]j, where k is odd . Thus n = mj +m 2 + . . . +m|j . 
With the summands in some fixed order, an edge two-coloring of 
the 2 n-gon, and therefore a ring graph, may be produced in the 
following manner. Start at some fixed edge in the 2n-gon and 
color the first m| edges black, proceeding clockwise from the 
starting edge. Color the next m 2 edges white, the next m 3 




(5.2) 



Given n obstacles, we begin with a partition of n into 
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edges black, etc., until all k summands have been used. The 
colors of the remaining n edges become fixed as this procedure 
is foil owed . 

Using the procedure given above, only partitions with 
an odd number of summands are used. This guarantees that the 
starting edge and the 2n^^ edge do not have the same color. 
Therefore, the two-colorings of the 2n-gon are not equivalent, 
and the associated ring'graphs are not isomorphic. Figure 5.7 
illustrates two equivalent two-colorings of an octagon (under 
rotation through 45 degrees) using two different partitions. 
The first partition of n uses an odd number of summands, the 
second partition uses an even number of summands. 



Start 




start 




2 , 



1 . 1 . 1 . 1 



Figure 5.7 Equivalent Two-colorings Resulting 
From Different Partitions 



It is possible that the number of candidate 
partitions, those with an odd number of summands, may actually 
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be less than N, the total number of non-isomorphic colorings 
calculated using Burnside's Lemma. When this occurs, there 
must exist at least one partition which, when ordered 
differently, is capable of producing more than one graph which 
is non-isomorphic to the others. Figure 5.8 shows, for n = 7, 
two non-equivalent two-colorings constructed from different 
orderings of the partition whose summands are 2, 2, 1,1,1. The 
question which arises'' is: How many non-equivalent two- 
colorings, and hence non-isomorphic ring graphs may be 
produced from a given partition of n into an odd number, k, of 
summands? 

Polya's Pattern Inventory Theorem may be used to 
answer this question. 



Start start 





Figure 5.8 Non-Equivalent Colorings Constructed 
Using the Summands 2, 2, 1,1,1 of n = 7 
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2. Polya's Pattern Inventory Theorem 

An inventory is a list which indicates how many items 
of a specific type exist [Ref. 7:p. 333]. This information 

may be encoded in the coefficients of a polynomial. Polya's 
Pattern Inventory Theorem follows. 

Let V be the set of vertices of a regular polygon. Let 
G be a group of permutations on the set V that acts to 
induce an equivalence relation on the colorings of V. The 
inventory of non-equivalent colorings of V using two 
colors is given by the generating function 

P(j ((b+w), (b^+w^),..., (b*' + w^)), 

where Pg is the cycle index for the group of symmetries of 
the polygon. The inventory using colors Cj^ , C 2 , . . . , c^ is 
given by 

mm m 

Pg{ £ C., £c/,..., £ C/) . [Ref. 7:p. 355] (5.3) 

j=l j=l j=l 

3. Establishing the Cycle Index 

We begin with the set of all ordered partitions of n 
into k summands, where k is an odd integer. Any two 
partitions will be considered to be equivalent if there exists 
a cyclic shift of the elements of one partition, a complete 
reversal of the elements of one partition, or a combination of 
these which makes the partitions identical. Thus the set of 
ordered partitions of n may be separated into equivalence 
classes, with each class producing a non-equivalent two- 
coloring of the 2n-gon. 

The cyclic shifts and reversals of elements in a 
partition of k summands correspond to the rotations of a k-gon 
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and reflections of a k-gon about an axis drawn through a 



vertex and bisecting its opposite edge. Therefore, it is 
reasonable to consider the number of non-equivalent vertex 
colorings of a k-gon where the number of colors used is equal 
to the number of different summands in the partition. For 
example, with n = 7, a partition of n might be 2, 2, 1,1,1. 
Thus, k = 5. Here the vertices of a pentagon would be colored 
using two colors since only two different summands (1 and 2) 
appear in the partition. 

In examining any k-gon where k is odd, it is seen that 
k rotations and k reflections are possible for a total of 2k 
rigid motions. The cycle index for the k-gon given these 
permutations is then produced. 

4. Applying Polya’s Theorem 

Once the cycle- index, Pq , has been established, 
Polya's Theorem may be applied, where the number of colors, m, 
is equal to the number of different summands in the given 
partition. After the cycle index has been expanded using 
these m colors, the coefficient of the term whose exponents 
are equal to the multiplicities of each of the integers in the 
partition of n is determined. This coefficient represents the 
number of non-isomorphic ring graphs which may be produced 
from the given partition of n. 

The sum of these coefficients, after applying this 
procedure to all partitions of n with an odd number of 
summands, is equal to the total number, N, of non-isomorphic 
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ring graphs given n fixed obstacles. This represents the 
number of distinct reference frames which may be constructed 
in the plane with n obstacles, thus generating N different 
sets of homotopy class names. 

E. AN ALTERNATE SOLUTION TO THE COUNTING PROBLEM 

Upon completing the research and derivation of the 
solution method discussed above, an alternate solution was 
encountered in the literature. 

Fredricksen and Kessler [Ref. 9] present a procedure which 



may be employed 


in counting 


and 


generating 


the set of non- 


isomorphic 


ring 


graphs . 


This 


method 


is 


based upon an 


algorithm 


which 


produces 


the 


set of 


al 


1 lexicographic 



compositions of a positive integer, n. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

Determination of the shortest path between two points in 
the plane containing obstacles has applications in many 
fields, particularly that of robotic path planning. This 
thesis presents a computational method which may be employed 
to generate the names of homotopy classes for a set of paths. 
By modeling the topological relationships of the region with 
a graph, a tool is made available for employing this method in 
producing only the names of those classes which contain 
candidates for the shortest path. 

B. RECOMMENDATIONS FOR FURTHER STUDY 

The field of graph theory is one that is gaining vast 
interest in mathematics and computer science. The graphs 
which are used to model the regions discussed in this thesis 
have several interesting characteristics. Further study and 
analysis of these graphs and their properties may aid future 
researchers in answering other questions concerning the 
shortest path problem. For instance, is there a sharp upper 
bound on the length of the class names which are generated by 
traversing the graph using the algorithm presented in 
Chapter III? 
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Another logical step in the shortest path problem is to 
extend the analyses presented here to include the third 
dimension. This extension follows since most industrial robot 
arms are of a pivoting arm type which are permitted to move 
freely through space while rotating at some fixed shoulder 
point . 

The computational investigation presented in Appendix A 
was programmed with knowledge of only elementary FORTRAN 
programming skills. Although a complexity analysis was not 
done, it is believed that further study and refinement of the 
program code could result in a program which is much more 
efficient . 

Finally, this thesis introduces only part of the solution 
to the shortest path problem. The thesis submitted by CPT 
Andre M. Cuerington, U.S. Army, presents the remaining 
portions of the solution, and may be studied in conjunction 
with this paper. 
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APPENDIX A. A COMPUTATIONAL INVESTIGATION 



A. INTRODUCTION 

This appendix describes a computational investigation 
which was conducted as an aid in verifying that the homotopy 
class names produced by Algorithm 1, are reasonable when 
compared to a known method of naming homotopy classes 
[Ref. 2]. 

A second algorithm. Algorithm 2, is introduced here to 
forge a link between the class names of Chapter II and the 
well-known fundamental group. Two processes are presented to 
exploit this link: If Algorithm 1 truly names the classes, 
then both procedures should produce the same output for every 
path tested. One million test cases were examined and no 
counter-examples were found. While this computational 
evidence is not a proof, it does support the proof of the 
claim that Algorithm 1 names homotopy classes [Ref. 2]. 

B. ALGORITHM 2 

1. Fundamental Group 

The class name obtained from Algorithm 1 was written 
in terms of the alphabet A defined in Chapter II. We present 
an algorithm here where the name obtained is expressed in 
terms of the fundamental group of a topological space T. 
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The basic idea for producing a fundamental group for 
this problem is to regard the paths in T as elements of the 
group, with path concatenation (*) as the group operation. 
There are two problems with considering the paths as the group 
elements. First, it is not necessarily possible to 
concatenate any two paths in T. In order to concatenate two 
paths using the operation (*), the first path must end at the 
point where the second begins. To make it possible to 
concatenate paths, in this analysis the test paths all begin 
and end at the point a. This point is called a base point. 

The second problem to overcome is that an inverse is 
not well defined for paths. If the identity is defined as the 
empty string, i.e., stay at the base point, then when a path 
and its inverse are concatenated we have traced out a path 
which is not equal to the identity. To avoid this problem, 
the homotopv classes of paths are considered. Then all paths 
that have an empty character string representation will be in 
the identity class. Also, any class followed by its inverse 
will equal the identity. In Figure A.l, gj represents the 
class / and (gj)*^ represents the class So, g^ * (gp'^ 
represents the class a2l3jB2a2 which reduces to the identity 
class by the rules of Algorithm 1. 

The group elements, therefore, are homotopy classes of 
loops around obstacles based at a common point, and the group 
operation is defined in terms of concatenation (*) of these 
classes . 
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Figure A.l A Homotopy Class and its Inverse 

Given a base point a, let p be any loop beginning and 
ending at a. Let [pj] denote the homotopy class of p^. So the 
set of group elements is {[pj] such that P' is a member of T, 
where pj is a loop based at a}. This set will be called G 
with the elements denoted gj. Figure A. 2 illustrates several 
paths and their homotopy classes given n = 2 obstacles. 




Figure A. 2 [pj] = [pj] * [P 3 ] 



Now the group operation (*) can be defined by 

[p] * [q] = [pq]- 
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It is important to note that the fundamental group is 



finitely generated. As Figure A. 3 illustrates, the generators 
are not unique. By fixing a set of generators, however, the 
class names become fixed. 



Figure A. 3 Alternate Generators For Fundamental Group of the 

Space With Two Obstacles 

With the above information describing the fundamental 
group representation. Algorithm 2 is given below. 




> 



a 
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2. Functions in Algorithm 2 

a. Side Array 

Consider the reference lines Lj^ constructed in 
Chapter II to be oriented rays in the direction from c to b|j . 
It then becomes reasonable to discuss a 'right' and a 'left' 
side of those lines. We allow a moving point to trace a path 
p from start point to destination point. Side(k) is a 
function which defines the side of Lj^ on which the moving 
point lies. The output is either 'left' or 'right'. The 
output is never 'on' because this routine is used only after 
a complete crossing of Lj^ takes place. A crossing is 
considered to be complete when the moving point leaves Ljj to 
one side after having met Lj^ from the opposite side. 

b. Switch Function 

The switch function is defined by: 

I 'right' if side (k) = 'left' 
switch (side(k)) = { 

I 'left' if side (k) = 'right' 

As the moving point traces p, each time a reference line is 

crossed, the switch function will be applied to indicate on 

which side of Lj^ the moving point lies. 

c. Index Function 

Let X|j be any character from the alphabet A 
representing the k^^ element in R(p) . Let index(x) = j if x 

equals Qj or Bj,’ thus, the k^^ crossing completed by the moving 
point is a crossing of Lj . So if the moving point is on the 



52 



left side of L 3 initially and this line is now crossed on 
either side of bj, then side(3) would equal 'right', 
d. The Algorithm 

Using the functions described above, Algorithm 2 
is given in Figure A. 4. The algorithm works in two phases. 
The first phase initializes the array side(k) for k = 1, . . . ,n. 
The second phase reads the raw string from left to right and 
adds an element to the fundamental group representation for 
each B crossing. When the end of the input string is reached 
the output, F(R(p)), is a shorter string with one character 
corresponding to every B element in the original raw string. 
Each character represents a generator or its inverse. 



begin 

input R(p) = xjXj . . . Xj, 
for j = 1 , . . . , n 

If a is to right of reference line Lj then 
side (j) - right 

else 

side (j) *- left 
end i f 
end for 
G - [a] 

for k = 1 , . . . , m 
i = index(X|j) 

side(i) - switch(side(i ) ) 
if X]^ = Bj for some r, then 

if side(index( X|j) ) = left, then 
G - G * g^ 

else 

G - G * (g:)‘l 
end if 
end i f 
end for 
F(R(p)) - G 
end 

Figure A. 4 Algorithm 2 
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e. Fundamental Group Cancellation Function 

The cancellation rules in the fundamental group 
differ somewhat from those for the raw string. Although sets 
of generators are not unique, we can obtain a unique 
representation of each class with respect to a particular set 
of generators. For any given set of generators, every 
homotopy class can be represented as a product of these 
generators and their inverses. Even so, this representation 
is not unique until cancellation is applied. The cancellation 
rule foil ows . 

Let G = {gp . . . , g^} be a set of generators of the 
fundamental group of T (base point a) and let Y = y^y 2 ---Y( be 
a representation of some homotopy class in terms of the gj and 
their inverses, i.e. for each i = 1 ,..., m, y^ = gj or {9)[)'^ 
for some j,k = 1 ,..., n. 

The cancellation function k is defined as follows. 
If Y contains a two character substring YiYiti with y^ = (y^^j)’^ 
then k(Y) is the string which results by removing both y,- and 
Y^^j, in their leftmost occurrence, from Y. Finitely many 
repeated applications of K produce a string in which no 
further cancellation is possible. We define this string as 
K(Y) . 
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C, THE COMPUTATIONAL INVESTIGATION 
1 . The Approach 

A graphical representation of the two procedures is 
provided in Figure A. 5. The investigation begins by accepting 
a random path p as input and generating two names for the 
homotopy class representing p which are output. 

The idea of the test is to determine for a given path 
its raw string character representation. For computational 
purposes, the paths considered are polygonal paths. This raw 
string is then input to both algorithms. 

The raw string (R(p)) is first entered into 
Algorithm 2. The output (F(R(p))) from Algorithm 2 is then 
subject to the cancellation routine which outputs the fully 
cancelled string (K(F(R(p) ) ) ) . Next, R(p) is input to 
Algorithm 1 which produces the canonical string (C(R(p))). 
Algorithm 2 is then applied to C(R(p)) to output F(C(R(p))). 
The results of these tests are then compared, and the program 
continues to run until a predetermined number of paths are 
checked . 

It can be shown that C(R(p)) is a unique class name if 
and only if F(C(R(p))) equals K(F(R(p))) for all p in 
T [Ref. 2]. Before the proof of this conjecture was obtained, 
the above procedures were programmed and tested for one 
million different paths and no counter-examples were 
discovered . 
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p 



ALG 2 



y 

R(p) 





k(F(R(p))) 



F(C(R(p))) 



K(F(R(p))) 



Figure A. 5 Flow Chart of Two Algorithms 



2 . The Test 

The general algorithm used to test the two procedures 
represented in Figure A. 5 is easily followed and is shown in 
Figure A. 6. However, the attached program, which was used to 
test the model, is substantially more involved. This was done 
in an attempt to write efficient programming code. 

Two time saving techniques are employed in order to 
test the one million cases in 25 minutes. The first method 
involves subroutines which perform cancellation. It is simple 
to code a program that scans the character representation 
repeatedly to find all possible cancellation in the string 

reduces to Bj. However, in the enclosed code. 
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pointers are used to mark the position where the first 
cancellation occurs and then check the newly adjacent 
characters for cancellation. This method reduces the 
complexity of an algorithm and is applied in both Algorithm 1 
and Algorithm 2. 



n = number of reference frames to be considered 

m = number of paths to be considered 

nobs = number of obstacles on each board 

seed = input seed for random generator (not needed?) 

begin 

for board = 1 , . . . , n 

create reference frame or board 
for path = 1 , . . . , n 

1. Create a polygonal path p with nseg 
segments 

2 . Form R( p ) 

3 . Form F(R(p) ) 

4. Form K(F(R(p) ) ) 

5 . Form C(R(p) ) 

6. Form F(C(R(p) ) ) 

if K(F(R(p))) does not equal F(C(R(p))) then 
print (raw string, board, and path info 
for the counter example) 
end i f 
end for 
end for 

print (final seed) 
end 



Figure A. 6 Test Algorithm 



The second time-saving device involves the sorting of 
character strings. Instead of a bubble sort method (order n 
complexity), a merge sort algorithm (order n(log(n)) 
complexity) is used. As can be seen in the attached code, 
this choice to save computer time calls for a significant 
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increase in programmer time, which is generally much more 
expensive . 

The code which is used to implement this test is 
presented in the remaining pages of this appendix. 
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PROGRAM CLASNAME FORTRAN 



C THIS PROGRAM RANDOM.Y GENERATES OBSTACLES AND A POLYGONAL PATH 
C THROUGH THOSE OBSTACLES TO DETERMINE WHETHER OR NOT THE PROCEDURES 
C USED BY COMPETING ALGOR ITFMS PRODUCE THE SAME FUNDAMENTAL GROUP 
C REPRESENTATION FOR THAT PATH. 

Q 

C THIS PORTION OF THE PROGRAM SERVES AS THE MAIN DRIVER WHICH RECEIVES 
C THE PARAMETERS DEFINING THE REGION AND INITIALIZES THE LINK LIST 
C ARRAYS WHICH WILL BE USED TO REPRESENT THE POLYGONAL PATH. 

INPUT: INITIAL SEED FOR THE RANDOM NUMBER GENERATOR, NUMBER OF 

OBSTACLE CONFIGURATIONS AND PATHS TO BE TESTED, NUMBER OF 
OBSTACLES IN THE PLANE AND NUMBER OF SEGMENTS IN EACH 
POLYGONAL PATH 

OUTPUT: FINAL RANDOM NUMBER GENERATOR SEED AND MESSAGES INDICATING 

ANY PATHS WHICH PRODUCE DIFFERENT FUNDAMENTAL GROUP 
REPRESENTATIONS 

REAL*8 BX(IOOO), BY(IOOO), X(1000), Y(1000), DSEED 
INTEGER NOBS, NSEQS, NUMPTS, N, M, HEAD(IOOO) 

+ NEXT(IOOO), PRED(IOOO), BOARD, PATH 

EXTERNAL GQUBFS 

DSEED = 123457.0 
M = 1000 
N = 1000 
NOBS = 20 
NSEQS = 5 
NUMPTS = NSEQS + 1 
PRINT*,’ INPUT SEED’, DSEED 

DO 1 BOARD = 1, N 

CALL BOARDS (DSEED, BX, BY, NOBS) 

DO 2 PATH = 1, M 

CALL IN IT (hEAD, NEXT, PRED, PSEED, DSEED ) 

CALL PATHS (DSEED, X, Y, NSEGS) 

CALL TEST (HEAD, NEXT, PRED, BX, BY, NOBS, X,Y, NUMPTS, BSEED, PSEED) 

2 CONTINUE 
1 CONTINUE 

PRINT*,’ ’ 

PRINT*, ’FINAL SEED’, DSEED 

STOP 
END 
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SUBROUTINE INIT ( HEAD, NEXT, PRED , PSEED, DSEED) 

C THIS SUBROUTINE INITIALIZES THE DOUBLE LINK LIST ARRAYS THAT WILL 
C REPRESENT THE PATH. 

C INPUT: DSEED 

C OUTPUT: HEAD, NEXT, AND PRED ARRAYS SET TO ZERO, AND SEED FOR THE 
C R.N.Q. PRIOR TO CONSTRUCTION OF THE PATH 

INTEGER h£AD(1000), NEXT(IOOO), PRED(IOOO) 

REAL*8 PSEED, DSEED 

PSEED = DSEED 

DO 3 K = 1,1000 
HEAD(K) = 0 
NEXT(K) = 0 
PRED(K) = 0 
3 CONTINUE 

RETURN 

END 



SUBROUTINE BOARDS (DSEED, BX, BY, NOBS) 

C THIS SUBROUTINE USES A PSEUDO RAhtfXDM NUMBER GENERATOR TO CREATE THE 
C COORDINATES OF EACH OBSTACLE ON THE BOARD, SCALING ALL COORDINATES 
C TO BE IN THE INTERVAL (-1,1). 

C INPUT: NUMBER OF OBSTACLES IN THE REGION AND A SEED FOR THE R.N.G. 

C OUTPUT: OBSTACLE COORDINATES 

REAL*8 BX(IOOO), BY(IOOO), DSEED 
EXTERNAL GGUBFS 

DO 1 I = 1, NOBS 

BX(I) = 2. * GQLBFS(DSEED) - 1. 

1 BY(I) = 2. * GQUBFS(DSEED) - 1. 

RETURN 

END 
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SUBROUTINE PATHS (DSEED, X, Y, NSEQS) 

Q llf^it******************************************************************* 

C RANDOMLY GENERATES THE X AND Y COORDINATES OF THE VERTICES FOR THE 
C POLYGONAL PATH, SCALING ALL COORDINATES TO BE IN TTC INTERVAL (-1,1). 
C THIS SUBROUTINE ALSO ENSURES THAT THE PATH IS A CLOSED LOOP BY 
C ASSIGNING THE START/FINISH POINTS THE SAME COORDINATES. 

C INPUT: NUMBER OF PATH SEGMENTS AND A SEED FOR THE R.N.G. 

C OUTPUT: COORDINATES OF VERTICES ALONG THE POLYGONAL PATH 

REAL*8 X(1000), Y(1000), DSEED 
EXTERNAL QGUBFS 

DO 1 1=1, NSEGS 

X(l) = 2. * QQUBFS(DSEED) - 1. 

1 Y(l) = 2. * Q3UBFS(DSEED) - 1. 

X(NSEQS+1) = X(1) 

Y(NSEQS+1) = Y(1) 

RETURN 

END 



SUBROUTINE TEST(HEAD, NEXT, PRED, BX, BY, NOBS, X, Y, I^JMPTS, 

+ BSEED, PSEED) 

Q 

C GENERATES THE RAW STRING OF CHARACTERS REPRESENTING A PATH MD 
C DETERMINES ITS FUNDA^ENTAL GROUP REPRESENTATION USING ALGORITHMS 1 
C AND 2. RESULTS OF THESE COMPETING ALGORITHMS ARE THEN COMPARED FOR 
C DIFFERENCES. 

C INPUT: A POLYGONAL PATH AND COORDINATES OF POINTS REPRESENTING 

C OBSTACLES 

C OUTPUT: COMPARED RESULTS OF ALGORITVMS 

INTEGER HEAD(IOOO), NEXT(IOOO), PRED(IOOO), FR(IOOO), NR(IOOO), 

+ PR(IOOO), FL(IOOO), NL(IOOO), PL(IOOO) 

REAL*8 BX(IOOO), BY(IOOO), X(1000), Y(1000) 

DATA FR/ 1000*0/ 

DATA NR/ 1000*0/ 

DATA PR/ 1000*0/ 

DATA FL/ 1000*0/ 

DATA NL/ 1000*0/ 

DATA PL/ 1000*0/ 
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GO GOO 



CALL RAWSTR(BX, BY, NOBS, X, Y, NUMPTS, HEAD, NEXT, PRED, NELEMS) 
CALL ALG2 (NOBS, NELEMS, HEAD, NEXT, X, Y, BX, BY, FR, NR, PR) 
CALL CANALQ2 (FR, NR, PR) 

CALL ALQ1 (HEAD, NEXT, PRED, NELEMS) 

CALL ALG2 (NOBS, NELEMS, HEAD, NEXT, X, Y, BX, BY, FL, NL, PL) 
CALL CTREX (FL, NL, PL, FR, NR, PR, BSEED, PSEED) 

RETURN 

END 



SUBROUTINE RAWSTR(BX, BY, NOBS, X, Y, NUMPTS, hCAD, NEXT, PRED, 

+ NELEMS) 

Q 

C PRODUCES THE RAW STRING OF CHARACTERS WHICH REPRESENTS THE PATH, 

C PAYING CLOSE ATTENTION TO THE ORDERING OF THE CHARACTERS WHERE 
C NECESSARY. THE STRING IS CONSTRUCTED BY IDENTIFYING THOSE OBSTACLE 
C REFERENCE LINES WHICH ARE CROSSED AS EACH PATH SEGMENT IS TRAVERSED 
C IN ORDER. 

INPUT: ORDERED LIST OF VERTICES REPRESENTING THE PATH AND A POINT 

REPRESENTING EACH OBSTACLE IN THE REGION 

OUTPUT: A RAW STRING OF ALPHAS AND BETAS CONTAINED IN AN ARRAY 

NAMED HEADO AND ITS PARALLEL ARRAYS NEXT() AND PRED() WHICH 
PRODUCE THE DOUBLE LINKED LIST 

IMPLICIT REAL*8 (A-H,0-Z) 

INTEGER HEAD, NEXT, SEG, FSTSTR, HDINDX, SEGEND, FIRST, START, 

+ LENGTH, PRED 

LOGICAL ALLALF 

DIMENSION X(1000), Y(1000), BX(IOOO), BY(IOOO), A1(1000), 

+ BI(IOOO), A2(1000), B2(1000), D2(1000), DIST(IOOO), 

+ FIRST(IOOO), HEAD(IOOO), NEXT(IOOO), PRED(IOOO) 

DATA D I ST/ 1000*0.0/ 

DATA FIRST/ 1000*0/ 

CALL SETUP(NUMPTS,NOBS,CX,CY,X,Y,BX,BY,A1,B1,A2,B2,D2) 

HDINDX = 2 
NSEGS = NUMPTS - 1 
LOLD = 1 
HEAD(LOLD) = 0 
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DO 7 SEQ = 1, NSEQS 
SEGEND = SEG + 1 
FSTSTR = HDINDX 
LSTSTR = HDINDX 
ALLALF = .TRUE. 
LENGTH = 0 



+ 



+ 



DO 6 LINE = 1, NOBS 

CHECK1 = Ui(LINE)*X(SEG)+B1(LINE)*Y(SEG))* 

(A1(LINE)*X(SEGEND)+B1(LINE)*Y(SEQEND)) 
IF (CHECK1.lt. 0) THEN 
LENGTH = LENGTH + 1 
CHECK2 = (A2(SEG)*BX(LINE)+B2(SEG)* 

BY( L I NE ) +D2 ( SEG) ) *D2( SEG) 

IF (CHECK2.lt. 0) THEN 
HEAD(HDINDX) = -LINE 
LSTSTR = HDINDX 
HDINDX = HDINDX + 1 



ELSE 

CALL CASES1 (A1 , B1 , A2, B2, D2, SEG, LINE, XINT, 
+ YINT) 

DISTC = XINT**2 + YINT«'*2 

DISTB = (XINT-BX(LINE))**2 + (YINT-BY(LINE))**2 
IF (DISTB.lt. DISTC) THEN 
HEAD(HDINDX) = LINE 
LSTSTR = HDINDX 
HDINDX = HDINDX + 1 
ALLALF = .FALSE. 

ELSE 

HEAD(HDINDX) = -LINE 
LSTSTR = HDINDX 
HDINDX = HDINDX + 1 
END IF 
END IF 
ENDIF 

6 CONTINUE 

IF (LENGTH. NE.O) THEN 
IF (ALLALF) THEN 

CALL ALPHAS (NEXT, LOLD, FSTSTR, LSTSTR, SEG, START, 

+ HEAD, NSEGS, PRED) 

ELSE 

CALL ORDER(SEG, LINE, HDINDX, HEAD, NEXT, A1 , A2, 

+ B1, B2, D2, FSTSTR, LSTSTR, X, Y, START, 

+ LENGTH, FIRST, LOLD, NSEQS, PRED) 

ENDIF 

ENDIF 

7 CONTINUE 



NELEMS) 



CALL COUNTR( START, NEXT, HEAD, 

RETURN 

END 
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SUBROUTINE SETUP(NUMPTS, NOBS, OX, CY, X, Y, BX, BY, A1, B1, A2, 

+ B2, D2) 

c 

C FOR EACH OBSTACLE, THIS ROUTINE DETERMINES THE COEFFICIENTS OF THE 
C EQUATION FOR THE REFERENCE LINE FROM THE OBSTACLE TO THE ORIGIN. IN 
C ADDITION, IT CALCULATES THE COEFFICIENTS OF THE LINE REPRESENTING 
C EACH SEGMENT OF THE POLYGONAL PATH. 

INPUT: NUMPTS, NOBS, COORDINATES OF VERTICES ALONG POLYGONAL PATH 

AND COORDINATES OF THOSE POINTS WHICH REPRESENT EACH 
OBSTACLE 

OUTPUT: COEFFICIENTS OF LINEAR EQUATIONS REPRESENTING PATH SEGMENTS 

AND REFERENCE LINES FOR EACH OBSTACLE 

IMPLICIT REAL*8 (A-H, O-Z) 

DIMENSION X(1000), Y(1000), BX(IOOO), BY(IOOO), A1(1000), 

+ BI(IOOO), A2(1000), B2(1000), D2(1000) 

DO 4 I = 1, NOBS 
Aid) = BY(I) 

Bid) = -BX(I) 

4 CONTINUE 

CX = 0.0 
CY = 0.0 
I = 1 

DO 5 J = 2, NUMPTS 

A2(l) = Y(J) - Y(l) 

B2(l) = X(l) - X(J) 

D2(l) = (Yd)*X(J))-(X(l)*Y(J)) 

I = 1 + 1 

5 CONTINUE 

RETURN 
END 
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SUBROLfTINE CASES1(A1, B1, A2, B2, D2, SEQ, LINE, XINT, YINT) 

C THIS SUBROUTINE DETERMINES THE COORDINATES FOR THE POINT OF 
C INTERSECTION OF A GIVEN PATH SEQhtNT AND A GIVEN OBSTACLE REFERENCE 
C LINE. NUMERICAL STABILITY OF CALCULATIONS REQUIRES THE MATHEMATICAL 
C OPERATIONS TO BE SEPARATED INTO CASES, THE BEST CASE BEING USED FOR 
C EACH PARTICULAR SITUATION. 

INPUT: COEFFICIENTS OF LINEAR EQUATIONS FOR PATH SEGMENT AND 

REFERENCE LINE TO BE EVALUATED 

OUTPUT: POINT OF INTERSECTION OF THE PATH SEQh€NT AND THE OBSTACLE 

REFERENCE LINE 

IMPLICIT REAL*8 (A-H,0-Z) 

INTEGER SEG 

DIMENSION A1(1000), B1(1000), A2(1000), B2(1000), D2(1000), A(2,2) 

A(l,l) = Al(LINE) 

AO. 2) = BI(LINE) 

A(2,1) = A2(SEG) 

A(2,2) = B2(SEG) 

BIGEST = 0.0 

DO 1 L = 1,2 

DO 1 K = 1,2 

TEST = DABS(A(K,L)) 

IF (TEST. GT. BIGEST) THEN 
BIGEST = TEST 
KBIG = K 
LBIG = L 
ENDIF 

1 CONTINUE 

IF (KBIG.EQ.1) THEN 

IF (LBIG.EQ.1) THEN 

YINT = -D2(SEG)/(B2(SEG)-B1(LINE)*A2(SEG)/A1(LINE)) 

XINT = -B1(LINE)*YINT/A1(LINE) 

RETURN 

ELSE 

XINT = -D2(SEG)/(A2(SEG)-A1(LINE)*B2(SEG)/B1(LINE)) 

YINT = -A1(LINE)*XINT/B1(LINE) 

RETURN 
ENDIF 

ELSE 

IF (LBIG.EQ.1) THEN 

YINT = (D2(SEG)*A1(LINE)/A2(SEG))/ 

+ (B1 (L I NE)-B2(SEG)*A1 (L I NE)/A2(SEG)) 

XINT = (-D2(SEQ)-B2(SEG)*YINT)/A2(SEG) 

RETURN 
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GO GO 



ELSE 



XINT = (D2(SEG)*B1(LINE)/B2(SEG))/ 

+ (A1(LINE)-A2(SEG)*B1(LINE)/B2(SEG)) 

YINT = (-D2(SEG)-XINT*A2(SEG))/B2(SEG) 
RETURN 
END IF 
ENDIF 

END 



SUBROUTINE ALPHAS (NEXT, LOLD, FSTSTR, LSTSTR, SEG, START, HEAD, 

+ NSEGS, PRED) 

Q t**********************ii^*****************t*********t*****t************ 

C GIVEN THAT A SEGMENT OF THE PATH CROSSES ONLY ALPHA RAYS WHEN 
C TRAVERSED, THESE CROSSINGS ARE SIMPLY INSERTED INTO THE LINK LIST 
C ARRAYS IN THE ORDER IN WHICH THEY WERE DETECTED, I.E., SMALLEST TO 
C LARGEST IN ABSOLUTE VALUE. 

INPUT; STRING REPRESENTING THE REFERENCE LINES CROSSED BY THE 

CURRENT PATH SEGMENT, GIVEN THAT ALL CROSSINGS ARE ALPHAS 

OUTPUT: UPDATED NEXT() AND PRED() ARRAYS WHICH CONTAIN THE STRING 

REPRESENTING THE MOST CURRENT PATH SEGMENT 

INTEGER NEXT (1000), FSTSTR, SEG, START, HEAD (1000), PRED (1000) 

NEXT(LOLD) = FSTSTR 
PRED (FSTSTR) = LOLD 
LAST = LSTSTR - 1 

IF ( LAST. GT. FSTSTR) THEN 
DO 12 I = FSTSTR, LAST 
NEXT(I) =1+1 
PRED(I+1) = I 
12 CONTINUE 

ELSE IF ( LAST. EQ. FSTSTR) THEN 
NEXT (FSTSTR) = LSTSTR 
PRED(LSTSTR) = FSTSTR 
ENDIF 

IF(SEG.EQ.I) START = LOLD 

LOLD = LSTSTR 
RETURN 
END 
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SUBROUTINE ORDER(SEG, LINE, HDINDX, HEAD, NEXT, Al, A2, B1, B2, 

+ D2, FSTSTR, LSTSTR, X, Y, START, LENGTH, FIRST, 

+ LOLD, NSEGS, PRED) 

C GIVEN THAT A SEGMENT OF THE PATH CROSSES ONE OR MORE BETAS, THIS 
C SUBROUTINE DETERMINES THE ACTUAL ORDER OF CROSSING WHEN TRAVERSING 
0 THE SEGMENT FROM BEGINNING TO END. THIS IS DONE BY FIRST DETERMINING 
C THE DISTANCE FROM THE SEGMENT START POINT TO THE POINT OF 
C INTERSECTION OF EACH CROSSED OBSTACLE REFERENCE LINE. THESE 
C DISTANCES ARE THEN SORTED FROM SMALLEST TO LARGEST AND CROSSINGS 
C ARE UPDATED IN THE LINK LIST ACCORDINGLY. 

C INPUT: HEAD() AND ALL COORDINATES REQUIRED TO DETERMINE THE 

C DISTANCES FROM INITIAL VERTEX OF PATH SEGMENT TO EACH OF 

C THE CROSSED REFERENCE LINES 

C OUTPUT: HEAD(), NEXT() AND PRED() ARRAYS CONTAINING THE RAW STRING 

C WHICH ACCURATELY REPRESENTS THE PATH ALONG THE CURRENT 

C SEGMENT 

INTEGER SEG, LINE, HEAD(IOOO), NEXT(IOOO), FSTSTR, HDINDX, 

+ START, FIRST(IOOO), PRED(IOOO), F 

REAL*8 XINTER(IOOO), YINTER(IOOO), A1(1000), Bl(IOOO), 

+ A2(1000), B2(1000), D2(1000), DIST(IOOO), X(1000), Y(1000) 

DO 8 J = FSTSTR, LSTSTR 
LINE = ABS(HEAD(J)) 

CALL CASES2 (Al, Bl, A2, B2, D2, SEG, LINE, XINTER, YINTER) 
DIST(J) = (XINTER(LINE)-X(SEG))**2+(YINTER(LINE)-Y(SEG))**2 
DIST(J) = -DIST(J) 

NEXT(J) = J+1 
PRED(J+1) = J 
8 CONTINUE 

NEXT (LSTSTR) = 0 
PRED(LSTSTR+1) = LSTSTR 

CALL MERG2(DIST, NEXT, PRED, FSTSTR, F) 

NEXT(LOLD) = F 
PRED(F) = LOLD 
IF(SEG.EQ.l) START = LOLD 
LOLD = NEXT (LOLD) 

17 IF(NEXT(LOLD).NE.O) THEN 
LOLD = NEXT (LOLD) 

GO TO 17 
ENDIF 

RETURN 

END 
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SUBROUTINE CASES2(A1 , B1 , A2, B2, D2, SEG, LINE, XINTER, YINTER) 

c 

C DETERMINES OOORDINATES FOR THE POINT OF INTERSECTION OF THE PATH 
C SEGMENT AND EACH OF THE REFERENCE LINES IT CROSSES. AGAIN, IN ORDER 
0 TO MAINTAIN NUMERICAL STABILITY, CALCULATIONS ARE MADE USING THE 
C MOST APPROPRIATE CLOSED FORM EQUATION. 

INPUT: COEFFICIENTS OF LINEAR EQUATIONS FOR CURRENT PATH SEGMENT 

AND ALL OBSTACLE REFERENCE LINES WHICH IT CROSSES 

OUTPUT: COORDINATES FOR POINTS OF INTERSECTION OF EACH OBSTACLE 

REFERENCE LINE WITH THE PATH SEGMENT 

ILLICIT REAL*8 (A-H, 0-2) 

INTEGER SEG 

DIMENSION A1(1000), B1(1000), A2(1000), B2(1000), D2(1000), 

+ XINTER(IOOO), YINTER(IOOO), A(2,2) 

A(1,1) = AI(LINE) 

A(1,2) = BI(LINE) 

A(2,1) = A2(SEG) 

A(2,2) = B2(SEG) 

BIGEST = 0.0 

DO 1 L = 1,2 

DO 1 K = 1,2 

TEST = DABS(A(K,L)) 

IF (TEST. GT. BIGEST) THEN 
BIGEST = TEST 
KBIG = K 
LBIG = L 
END IF 

1 CONTINUE 

IF (KBIG.EQ.1) THEN 

IF (LBIG.EQ.1) THEN 

YINTER(LINE)= -D2(SEG)/(B2(SEG)-B1(LINE)*A2(SEG)/A1(LINE)) 
XINTER(LINE)= -B1(LINE)*YINTER(LINE)/A1(LINE) 

RETURN 

ELSE 

XINTER(LINE)= -D2(SEG)/(A2(SEG)-A1 (L INE)*B2(SEG)/B1 (L I NE) ) 
YINTER(LINE)= -A1 (L I NE)*XI NTER(L I NE)/B1 (L INE) 

RETURN 
ENDIF 

ELSE 

IF (LBIG.EQ.1) THEN 

YINTER(LINE) = (D2(SEG)*A1 (L I NE)/A2(SEG) )/ 

+ (B1(LINE)-B2(SEG)*A1(LINE)/A2(SEG)) 
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XINTER(LINE) = (-D2(SEG)-B2(SEG)*YINTER(L)NE))/A2(SEG) 
RETURN 

ELSE 

XINTER(LINE) = (D2(SEG)*B1 (L I NE)/B2(SEG) )/ 

+ (A1(LINE)-A2(SEG)*B1(LINE)/B2(SEG)) 

YINTER(LINE) = (-D2(SEG)-XINTER(LINE)*A2(SEG))/B2(SEG) 
RETURN 
END IF 
EfOIF 
END 



SUBROUTINE MERG2 (DIST, NEXT, PRED, FSTSTR, F) 

0 4;4;:t:*4:3f:****4:**4:*******4:3t;****4:***4:*t**:|c*4;4c:|c:tc4c***t4:4:***********4:4:****#**3f: 

C THIS SUBROUTINE SORTS A SUBSTRING OF ALL POSTIVE INTEGERS INTO 
C INCREASING ORDER. 

C INPUT: A DOUBLE LINK LIST CONSISTING OF DIST, NEXT, AND PRED ARRAYS 

OUTPUT: A DOUBLE LINK LIST WITH ALL ENTRIES PLACED IN THE ORDER 
IN WHICH THEIR RESPECTIVE REFERENCE LINES WERE CROSSED 

IMPLICIT INTEGER(A-Z) 

REAL*8 DIST(IOOO) 

LOGICAL DONE 

DIMENSION NEXT(IOOO), PRED(IOOO) 

DONE = .FALSE. 

P = FSTSTR 

1 F = P 

CALL SORT2 (F, PREDF, P, DIST, NEXT, DONE, PRED) 

IF(DONE) RETURN 
GOTO 1 
END 
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SUBROLfTINE SORT2(F, PREDF, P, DIST, NEXT, DONE, PRED) 

C MARCHES DOWN A SUBSTRING OF POSITIVE INTEGERS AND ONLY SORTS IF 
C ELEMENTS ARE NOT IN INCREASING ORDER. IF A NUMBER NEEDS TO BE 
C PLACED HIGHER IN THE LIST SUBROUTINE 'PUT' IS CALLED TO DO SO 

INPUT: POINTER F INTO DOUBLE LINK LIST ARRAYS DIST, NEXT, AND PRED 

TO INDICATE THE BEGINNING OF A SUBSTING OF POSITIVE INTEGERS, 
AND PREDF 

OUTPUT: P IS A POINTER, DIST(P) IS THE LAST POSITIVE INTEGER IN THE 
SUBSTRING THAT IS BEGUN BY DIST(F) 

IMPLICIT INTEGER (A-Z) 

REAL*8 DIST (1000) 

LOGICAL DONE 

DIMENSION NEXT(IOOO), PRED(IOOO) 

TAIL = F 

1 NTAIL = NEXT(TAIL) 

IF(NTAIL.EQ.O) THEN 
DONE = .TRUE. 

RETURN 
END IF 

IF(DIST(NTAIL).LE.DIST(TAIL)) THEN 
TAIL = NTAIL 
GOTO 1 
END IF 

CALL PUT2(F, TAIL, NTAIL, PREDF, DIST, NEXT, PRED) 

GOTO 1 
END 
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SUBROUTINE PUT2 (F, TAIL, NTAIL, PREDF, DIST, NEXT, PRED) 



C REARRANGES POINTERS TO PLACE POSITIVE INTEGERS IN INCREASING ORDER. 



C INPUT: F 

C TAIL 

C NTAIL 

C 
C 

C PREDF 



START OF POSITIVE SUBSTRING 
END OF CURRENTLY SORTED PORTION OF SUBSTRING 
NEXT(TAIL) - POINTER TO THE SUCCESSOR OF TAIL IN THE 
ITEM TO BE INCORPORATED INTO THE SORTED PORTION OF 
THE LIST 

THE PREDECESSOR OF F IN LINKED LIST 



OUTPUT: SOME POINTERS IN NEXT AND PRED ARE CHANGED TO PUT DIST(NTAIL) 
IN ITS PROPER PLACE IN THE SORTED PORTION ON THE LIST 



IMPLICIT INTEGER (A-2) 

REALMS DIST (1000) 

DIMENSION NEXT(IOOO), PRED(IOOO) 

IF (DIST(NTAIL).GE.DIST(F)) THEN 
TEMP = NEXT(NTAIL) 

NEXT(NTAIL) = F 
PRED(F) = NTAIL 
NEXT(TAIL) = TEMP 
PRED (TEMP) = TAIL 
F = NTAIL 
RETURN 
END IF 

CALL WALK2(F, TAIL, NTAIL, DIST, NEXT, PRED) 

RETURN 

END 
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SUBROUTINE WALK2(F, TAIL, NTAIL, DIST, NEXT, PRED) 

C 

THIS SUBROUTINE WALKS DOWN THE LINKED LIST AND PLAOES THE DIST(NTAIL) 
IN THE CORRECT POSITION IN THE DIST ARRAY. 

INPUT: F, TAIL, NTAIL ARE USED AS ABOVE 

OUTPUT: ALTERS POINTERS IN NEXT AND PRED ARRAYS TO PLACE DIST (NTAIL) 
AFTER DIST(F) AND BEFORE DIST(TAIL) IN DOUBLE LINK LIST 

IMPLICIT INTEGER (A-Z) 

REAL*8 DIST (1000), HNTAIL 
DIMENSION NEXT(IOOO), PRED(IOOO) 

I = F 

NEXTI = NEXT(I) 

HNTAIL = DIST(NTAIL) 

1 IF(HNTAIL.QE.DIST(NEXTI )) THEN 
NEXT(I) = NTAIL 
PRED(NTAIL) = I 
NNTAIL = NEXT(NTAIL) 

NEXT(NTAIL) = NEXTI 
PRED(NEXTI) = NTAIL 
NEXT(TAIL) = NNTAIL 
PRED(NNTAIL) = TAIL 
RETURN 
END IF 

I = NEXTI 
NEXTI = NEXT(I) 

GOTO 1 

END 
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SUBROUTINE OOUNTR (START, NEXT, HEAD, NELEHS) 

C THIS SUBROUTINE COUNTS THE NUMBER OF ELEMENTS IN A GIVEN STRING OF 
C CHARACTERS 

C INPUT: HEADO AND NEXT() ARRAYS FOR THE STRING OF CHARACTERS 

C OUTPUT: NUMBER OF ELEMENTS IN THE STRING 

INTEGER PTR, NEXT(IOOO), HEAD(IOOO), START 

PTR = NEXT(START) 

NELEMS = 0 

20 IF (PTR.NE.O) THEN 

NELEMS = NELEMS + 1 
PTR = NEXT (PTR) 

GO TO 20 
END IF 

RETURN 

END 
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SUBROUTINE ALQ2(NOBS, NELEMS, HEAD, NEXT, X, Y, BX, BY, FR, NR, 
+ PR) 

C GENERATES THE WELL-KNOW FUNDAMENTAL GROUP REPRESENTATION OF AN 
C EQUIVALENCE CLASS. 

INPUT: RAW OR CANONICAL STRING, WITH # OBSTACLES, # ELEMENTS IN 

STRING, COORDINATES OF A AND OBSTACLES (BK) 

OUTPUT: 1) IF INPUT IS RAW STRING, THEN POSSIBLY UNREDUCED 

FUNDAMENTAL GROUP WILL RESULT (IF CANCELLATION OF LIKE 
POSITIVE NUMBERS COULD HAVE OCCURRED IN THE RAW STRING). 

2) IF INPUT IS CANONICAL STRING, 1TCN THE RESULTING 
FUNDAMENTAL GROUP WILL BE IN REDUCED FORM. 

REAL*8 XA, YA, XB, YB, X(1000), Y(1000), BX(IOOO), BY(IOOO) 
INTEGER NOBS, NELEMS, HEAD(IOOO), NEXT(IOOO), PRED(IOOO), 

+ FR(IOOO), NR(IOOO), PR(IOOO), S(1000), FUNDGP ( 1 000 ) , 

+ START 

LOGICAL RIQHT(IOOO), RITE 

LENGTH = 1 
START = NEXT(1) 

XA = X(1) 

YA = Y(1) 

DO 1 K = 1, NOBS 
XB = BX(K) 

YB = BY(K) 

1 RIGHT(K) = RITE(XA, YA, XB, YB) 

NELEM = NELEMS + 1 
DO 10 M = 2, NELEM 

S(M) = HEAD(START) 

J = IABS(S(M)) 

RIGHT(J) = .NOT. RIGHT(J) 

IF(S(M).GT.O) THEN 

LENGTH = LENGTH + 1 
IF(RIGHT(J)) THEN 

FUNDGP (LENGTH) = J 

ELSE 

FUNDGP (LENGTH) = -J 

END IF 
ENDIF 

START = NEXT( START) 

IF(START.EQ.O) GOTO 11 
10 CONTINUE 
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11 FR(1) = 0 
NR(1) = 2 

DO 2 I = 2, LENGTH 
FR(I) = FUNDQP(I) 
NR(I) = I + 1 
PR(I) = I - 1 
2 CONTINUE 

NR(LENGTH) = 0 

RETURN 

END 



LOGICAL FUNCTION RITE(XA, YA, XB, YB) 

C THIS FUNCTION DETERMINES WHICH SIDE OF A GIVEN DIRECTED LINE ANY 
C POINT LIES. 

C INPUT: TWO POINTS THAT DETERMINE THE LINE 

OUTPUT: LOGICAL VARIABLE THAT IS TRUE IF A POINT LIES TO THE RIGHT 
AND FALSE IF A POINT LIES TO THE LEFT 

REAL*8 XB, YB, XA, YA, SIGNA 

RITE = .TRUE. 

SIGNA = -((YB*XA)-(XB*YA)) 

IF(SIGNA.LT.O)THEN 
RITE = .FALSE. 

RETURN 
END IF 

IF(SIGNA.EQ.O) THEN 

PRINT*, 'THE POINT A LIES ON THE LINE LK PROGRAM STOPS' 

STOP 
END IF 

RETURN 
END 
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SUBROUTINE CANALQ2(HEAD, NEXT, PRED) 

C THIS SUBROUTINE TAKES AN UNREDUCED FUNDAMENTAL GROUP REPRESENTATION 
C OF A GIVEN CLASS AND CANCELS A GENERATOR IF IT IS ADJACENT TO ITS 
C INVERSE. 

C INPUT: HEAD, NEXT, AND PRED ARRAYS 

OUTPUT: HEAD, NEXT, AND PRED ARRAYS WITH NEXT AND PRED REARRANGED 
TO SKIP AROUND CANCELLED ELEMENTS 

IMPLICIT INTEGER (A-Z) 

DIMENSION HEAD(IOOO), NEXT(IOOO), PRED(IOOO) 

START = 1 
PTR1 = START 
PTR2 = START 
PTR3 = NEXT(PTR2) 

10 IF (PTR3.NE.0) THEN 

IF (HEAD(PTR2).EQ.-(HEAD(PTR3))) THEN 
NEXT(PTRI) = NEXT(PTR3) 

PRED(NEXT(PTR3)) = PTR1 
IF(NEXT(PTR1).EQ.O) RETURN 
PTR2 = PTR1 
PTRl = PRED(PTRI) 

PTR3 = NEXT(PTR2) 

ELSE 

PTRl = PTR2 
PTR2 = PTR3 
PTR3 = NEXT(PTR2) 

END IF 
GO TO 10 
END IF 
RETURN 
END 
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SUBROUTINE ALG1(HEAD, NEXT, PRED, NELEMS) 

Q 

C THIS SUBROUTINE TAKES A GIVEN RAW STRING OF CHARACTERS REPRESENTING 
C A PATH AND PRODUCES THE CANONICAL FORM OF THAT STRING. THIS IS DONE 
C BY FIRST ORDERING ALL OF THE ALPHA SUBSTRINGS FROM SMALLEST TO 
C LARGEST IN ABSOLUTE VALUE. NEXT, CANCELLATION IS PERFORMED TO 
C ELIMINATE ALL LIKE PAIRS OF ADJACENT ELEMENTS FROM THE STRING. 

INPUT: RAW STRING IN FORM OF DOUBLE LINKED LIST WITH HEAD(), NEXT() 

AND PREDO arrays 

C OUTPUT: CANONICAL FORM OF THE RAW STRING 

ir^LICIT INTEGER (A-Z) 

DIMENSION HEAD(IOOO), NEXT(IOOO), PRED(IOOO) 

START = 1 

CALL MERG1(HEAD, NEXT, PRED) 

CALL CANCEL (START, HEAD, NEXT, PRED, NELEMS) 

RETURN 

END 
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SUBROUTINE MERG1 (HEAD, NEXT, PRED) 

Q ****:(c:)c:t:******:(c******************************************************** 

C THIS SUBROUTINE DOES THE INITIAL SORT OF THE ALPHA SUBSTRINGS IN THE 
C RAWSTRING (INCREASING IN ABSULUTE VALUE). 

INPUT: HEAD, NEXT, PRED ARRAYS REPRESENTING A DOUBLE LINKED LIST 

OF THE RAWSTRING WITH THE ALPHA SUBSTRINGS UNORDERED 

0 OUTPUT: POINTERS STORED IN THE ARRAYS NEXT AND PRED ARE ALTERED SO 
C THAT EACH SUBSTRING OF THE STORED LIST WHICH CONSISTS 

C ENTIRELY OF NEGATIVE INTEGERS IS SORTED INTO NON- I NCREASI NG 

C ORDER, WHILE SUBSTRINGS OF POSITIVE INTEGERS ARE LEFT 

C UNALTERED. 

ir^LICIT INTEGER(A-Z) 

LOGICAL DONE 

DIMENSION HEAD(IOOO), NEXT(IOOO), PRED(IOOO) 

DONE = .FALSE. 

P = 1 

1 CALL FRONT (P, F, PREDF, HEAD, NEXT, DONE) 

IF(DONE) RETURN 

CALL SORT (F, PREDF, P, HEAD, NEXT, DONE, PRED) 

IF(DONE) RETURN 

GOTO 1 

END 
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SUBROOTINE FRONT(P, F, PREDF, HEAD, NEXT, DONE) 

C FINDS THE BEGINNING OF NEGATIVE INTEGER STRINGS (ALPHA STRING) 

C INPUT: POINTER P INTO LINKED LIST 

OUTPUT: F IS A POINTER INTO THE HEAD ARRAY. POINTS TO FIRST NEGATIVE 
ENTRY WHICH OCCURS STRICTLY AFTER HEAD(P). PREDF IS POINTER 
SUCH THAT HEAD(PREDF) IS THE PREDESSOR OF HEAD(F). 

IMPLICIT INTEGER(A-Z) 

LOGICAL DONE 

DIMENSION HEAD(IOOO), NEXT(IOOO) 

F = P 

1 PREDF = F 
F = NEXT(F) 

IF (F.EQ.O) THEN 
DONE = .TRUE. 

RETURN 
ENDIF 

IF (HEAD(F).LT.O) RETURN 
GOTO 1 
END 
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SUBROUTINE SORT(F, PREDF, P, HEAD, NEXT, DONE, PRED) 

C MARCHES DOWN A SUBSTRING OF NEGATIVE INTEGERS AND ONLY SORTS IF 
C ELEMENTS ARE IN INCREASING ORDER. IF A NUh«ER NEEDS TO BE 
C PLACED HIGHER IN THE LIST SUBROLfTINE ’PUT’ IS CALLED TO DO SO 

INPUT: POINTER F INTO DOLBLE LINK LIST ARRAYS HEAD, NEXT, AND PRED 

TO INDICATE THE BEGINNING OF A SUBSTRING OF NEGATIVE 
INTEGERS; AND PREDF 

OinPUT: PISA POINTER, HEAD(P) IS THE LAST NEGATIVE INTEGER IN THE 
SUBSTRING THAT IS BEGUN BY HEAD(F) 

IMPLICIT INTEGER (A-Z) 

LOGICAL DONE 

DIMENSION HEAD(IOOO), NEXT(IOOO), PRED(IOOO) 

TAIL = F 

1 NTAIL = NEXT(TAIL) 

IF(NTAIL.EQ.O) THEN 
DONE = .TRUE. 

RETURN 
END IF 

IF(HEAD(NTAIL).GT.O) THEN 
P = TAIL 
RETURN 
END IF 

IF(HEAD(NTAIL).LE.HEAD(TAIL)) THEN 
TAIL = NTAIL 
GOTO 1 
END IF 

CALL PUT(F, TAIL, NTAIL, PREDF, HEAD, NEXT, PRED) 

GOTO 1 
END 
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SUBROUTINE PUT (F, TAIL, NTAIL, PREDF, HEAD, NEXT, PRED) 

C SUBROUTINE THAT REARRANGES POINTERS TO PLACE A NEGATIVE INTEGER IN 
C NON- INCREASING ORDER. 



0 

C 

C 

C 

C 

C 



INPUT: 



F - START OF NEGATIVE SUBSTRING 
TAIL - END OF CURRENTLY SORTED PORTION OF SUBSTRING 
NTAIL = NEXT(TAIL) - POINTER TO THE SUCCESSOR OF TAIL IN THE 
ITEM TO BE INCORPORATED INTO THE SORTED PORTION OF 
THE LIST. 

PREDF - THE PREDECESSOR OF F IN LINKED LIST 



C OUTPUT: SOME POINTERS IN NEXT AND PRED ARE CHANGED TO PUT HEAD(NTAIL) 
C IN ITS PROPER PLACE IN THE SORTED PORTION ON THE LIST 



IMPLICIT INTEGER (A-Z) 

DIMENSION HEAD(IOOO), NEXT(IOOO), PRED(IOOO) 

IF (HEAD(NTAIL).GE.HEAD(F)) THEN 
NEXT (PREDF) = NTAIL 
PRED (NTAIL) = PREDF 
TEMP = NEXT( NTAIL) 

NEXT(NTAIL) = F 
PRED(F) = NTAIL 
NEXT(TAIL) = TEMP 
PRED(TEMP) = TAIL 
F = NTAIL 
RETURN 

ENDIF 

CALL WALK(F, TAIL, NTAIL, HEAD, NEXT, PRED) 

RETURN 

END 
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SUBROUTINE WALK(F, TAIL, NTAIL, HEAD, NEXT, PRED) 

Q tttt*******^******************************************************^*** 

C THIS SUBROUITNE WALKS DOWN THE LINKED LIST AND PLACES THE HEAD(NTAIL) 

C IN THE CORRECT POSITION IN THE HEAD ARRAY. (DECREASING ORDER) 

C INPUT: F, TAIL, NTAIL ARE USED AS ABOVE 

C OUTPUT: ALTERS POINTERS IN NEXT AND PRED ARRAYS TO PLACE HEAD( NTAIL) 

C AFTER HEAD(F) AND BEFORE HEAD(TAIL) IN DOUBLE LINK LIST 

IMPLICIT INTEGER (A-Z) 

DIMENSION HEAD(IOOO), NEXT(IOOO), PRED(IOOO) 

I = F 

NEXT I = NEXT(I) 

HNTAIL = HEAD(NTAIL) 

1 IF(HNTAIL.GE.HEAD(NEXTI )) THEN 
NEXT(I) = NTAIL 
PRED(NTAIL) = I 
NNTAIL = NEXT(NTAIL) 

NEXT(NTAIL) = NEXT I 
PRED (NEXT I ) = NTAIL 
NEXT(TAIL) = NNTAIL 
PRED( NNTAIL) = TAIL 
RETURN 
END IF 

I = NEXT I 
NEXT I = NEXT(I) 

GOTO 1 

END 



82 



SUBRCHJTINE CANCEL (START, HEAD, NEXT, PRED, NELEMS) 

Q *******************************************************:^************:4e* 

C GIVEN THE SORTED RAW STRING, THIS SUBROUTINE REDUCES THE STRING TO 
C CANONICAL FORM BY ADJUSTING THE LINK LIST TO SKIP ANY PAIRS OF 
C ADJACENT LIKE ELEMENTS IN THE STRING. ONCE TWO ELEMENTS ARE 
C ELIMINATED FROM THE STRING, THE SUBSTRINGS WHICH WERE ON THEIR LEFT 
C AND RIGHT CONCATENATE TO FORM A NEW STRING AND HENCE A NEW PAIR OF 
C ADJACENT ELEMENTS WHICH MUST ALSO BE CHECKED FOR EQUALITY. IN THE 
C EVENT A BETA SUBSTRING BECOMES ANNHILATED FROM THE STRING, THE TWO 
C ADJACENT ALPHA STRINGS CONCATENATE AND ARE AGAIN SORTED FROM SMALLEST 
C TO LARGEST IN ABSOLUTE VALUE AS A SINGLE SUBSTRING. 

C INPUT: SORTED RAW STRING 

C OUTPUT: CANONICAL FORM OF THE RAW STRING 

IMPLICIT INTEGER (A-Z) 

DIMENSION HEAD(IOOO), NEXT(IOOO), PRED(IOOO) 

START = 1 
START 1 = START 
START2 = START 
BETA = START 
PTR1 = START 1 
PTR2 = START 1 
50 PTR3 = NEXT(PTR2) 

IF (HEAD(PTR3).LT.O) THEN 
START 1 = PTR3 
START2 = START 1 

ELSE 

BETA = PTR3 
END IF 

10 IF (PTR3.NE.0) THEN 

CALL CHEKER (HEAD, START1, START2, PTR1, PTR2, PTR3, BETA) 

IF (HEAD(PTR2).EQ.HEAD(PTR3)) THEN 
NEXT(PTRI) = NEXT(PTR3) 

PRED(NEXT(PTR3)) = PTR1 
PTR2 = NEXT(PTR3) 

IF (PTR2.EQ.0) GO TO 60 
CHECK1 = HEAD(PTR2) * HEAD(PTR3) 

CHECK2 = HEAD(PTRI) * HEAD(PTR2) 

IF (CHECK 1.LT.0. AND. CHECK2.GT.0. AND. HEAD(PTR2).LT.0)THEN 
PTR2 = PTR1 
PTR1 = PRED(PTRI) 

PTR3 = NEXT(PTR2) 

START 1 = START2 

6 IF(START1.NE. 1 .AND.HEAD(PRED(START1)) .LT.0)THEN 
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START 1 = PRED( START 1) 

00 TO 6 
ENDIF 

START2 = PTR3 

CALL MERGE (HEAD, NEXT, START1, START2, PRED) 

PTR1 = PRED ( START 1) 

PTR2 = PTR1 
00 TO 50 

ELSE IF (CHECK 1.LT.0. AND. CHECK2.EQ.0. AND. HEAD(PTR2).LT.O) 
+ THEN 

PTR2 = PTR1 
PTR3 = NEXT(PTRI) 

START1 = PTR3 
START2 = START 1 
GO TO 10 

ELSEIF (CHECK1 .GT.0.AND.CHECK2.GT.0.AND.HEAD(PTR2) .GT.O) 
+ THEN 

IF (BETA.EQ.2) THEN 
BETA = 1 

ELSEIF (BETA.EQ.1) THEN 
GO TO 9 

ELSE 

IF(STARTI.NE.I) BETA = PRED(STARTI) 

11 IF(HEAD(PRED(BETA)).GT.O) THEN 

BETA = PRED (BETA) 

GO TO 11 
ENDIF 
ENDIF 

9 PTR1 = PRED (BETA) 

PTR2 = BETA 
PTR3 = NEXT(PTR2) 

GO TO 10 

ELSEIF (CHECK1.LT.0. AND. CHECK2. GT.O. AND. HEAD(PTR2). GT.O) 
+ THEN 

PTR2 = PTR1 
PTR1 = PRED(PTRI) 

IF(HEAD(PTR1).LT.0)THEN 
START 1 = PTR1 

5 IF (HEAD(PRED(START1)).LT.0)THEN 

START1 = PRED(STARTI) 

GO TO 5 

ENDIF 

ENDIF 

PTR3 = NEXT(PTR2) 

GO TO 10 

ELSEIF (CHECK1.LT.O.AND.CHECK2.LT.O.AND.HEAD(PTR2).LT.O) 
+ THEN 

IF(STARTI.EQ.I) THEN 
START 1 = PTR2 
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START2 = PTR2 
PTR3 = NEXT(PTR2) 

ELSE 

START 1 = START2 
PTR3 = NEXT(PTR2) 

END IF 
00 TO 10 

ELSE 

START 1 = START2 
PTR3 = NEXT(PTR2) 

END IF 

ELSE 

PTR1 = PTR2 
PTR2 = PTR3 
PTR3 = NEXT(PTR2) 

END IF 
00 TO 10 
END IF 

60 CALL COUNTR (START, NEXT, HEAD, NELEMS) 

RETURN 

END 
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SUBROUTINE CHEKER(HEAD, START 1 , START2, PTR1, PTR2, PTR3, BETA) 

C 

C AS POINTERS MOVE ALONG THE CHARACTER STRING DURING CANCELLATION TMIS 
C SUBROUTINE CHECKS TO DETERMINE WHETHER OR NOT THE END OF ONE 
C SUBSTRING IS REACHED AND A NEW ONE BEGINS. DEPENDING ON THE OUTCOME 
C OF THIS CHECK, THE POINTERS USED TO IDENTIFY THE BEGINNING OF THE 
C TWO LATEST ALPHA STRINGS AND THE LATEST BETA STRING ARE UPDATED 

INPUT: HEAD() ARRAY AND PRESENT POINTER LOCATIONS FROM 'CANCEL' 

SUBROUTINE 

C OUTPUT: MODIFIED INDICES FOR LOCATION OF ALPHA AND BETA STRINGS 

INTEGER HEAD(IOOO), START 1 , START2, PTR1 , PTR2, PTR3, CHECK, BETA 
CHECK = HEAD(PTR3) ♦ HEAD(PTR2) 

IF (CHECK. GT.O) THEN 

IF (HEAD(PTR2). LT.O. AND. HEAD(PTRI). GT.O) START2 = PTR2 
IF (HEAD(PTR2). GT.O. AND. HEAD(PTRI). LT.O) BETA = PTR2 
ELSE IF (CHECK. LE.O) THEN 

IF (HEAD(PTR2). GT.O. AND. HEAD(PTRI). LT.O) THEN 
BETA = PTR2 
START 1 = START2 
START2 = PTR3 

ELSE IF (HEAD(PTR2).GT.O) THEN 
IF (START1.NE. 1) THEN 
START 1 = START2 
START2 = PTR3 

ELSE 

START 1 = PTR3 
START2 = START 1 
END IF 

ELSEIF (HEAD(PTR2). LT.O. AND. HEAD(PTRI). GT.O) THEN 
BETA = PTR3 
START 1 = START2 
START2 = PTR2 

ELSEIF (HEAD(PTR2). LT.O. AND. HEAD(PTRI). LT.O) THEN 
BETA = PTR3 
END IF 
END IF 

RETURN 

END 
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SUBROUTINE MERGE (HEAD, NEXT, ST ART 1, START2, PRED) 

C THIS SUBROUTINE MERGES TWO ORDERED ALPHA STRINGS. THE RESULTING 
C SUBSTRING WILL BE INCREASING IN ABSOLUTE VALUE. 

INPUT: HEADO, NEXT() AND PRED() ARRAYS ALONG WITH INDICES START 1 

AND START2 INDICATING BEGINNING OF THE ALPHA STRINGS TO BE 
MERGED 

OUTPUT: MODIFIED STRING CONTAINING THE MERGED ALPHA STRINGS WHICH 

WERE A RESULT OF ANNHILATION OF A BETA STRING 

INTEGER HEAD(IOOO), NEXT(IOOO), START1, START2, PI, P2, TAIL, 

+ PRED(IOOO) 

LOGICAL GO 

IF (START1.EQ.1) GO TO 2 

IF(HEAD(PRED(START1)).EQ.O) START1 = NEXT(PRED(START1 ) ) 

I F ( NEXT( PRED( START1 ) ) . NE . START 1 )THEN 
START 1 = NEXT(PRED(START1)) 

7 IF (START1.EQ. 1) GO TO 2 
I F (HEAD( START 1 ) . GT . 0 )THEN 

START 1 = NEXT ( START 1) 

GO TO 7 
END IF 

8 IF(HEAD(PRED(START1)).LT.0)THEN 
START 1 = PRED( START 1) 

GO TO 8 

END IF 
END IF 

2 CALL SET(HEAD, NEXT, START1, START2, PI, P2, TAIL, PRED) 

1 IF(GO(P1, P2, START2, HEAD, NEXT)) THEN 

CALL MERGER(HEAD, NEXT, PI, P2, TAIL, PRED) 

GOTO 1 
END IF 

CALL FXTAIL(HEAD, NEXT, PI, P2, TAIL, START2, PRED) 

RETURN 
END 
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SUBROUTINE SET(HEAD, NEXT, START 1 , START2, PI, P2, TAIL, PRED) 



C SETS THE STARTING POINTER OF THE NEW STRING AND THE TWO POINTERS 
C OF THE STRINGS TO BE MERGED 

C INPUT: POINTERS INTO TWO ALPHA STRINGS 

C OUTPUT: POINTER INTO BEGINNING OF A NEW SORTED STRING 

INTEGER HEAD(IOOO), NEXT(IOOO), START1, START2, PI, P2,TAIL, 

+ PRED(IOOO) 

PI = START 1 
P2 = START2 

IF(HEAD(P1).GT.HEAD(P2)) THEN 
TAIL = PI 
START 1 = PI 
PI = NEXT(PI) 

ELSE 

TAIL = P2 
START 1 = P2 
NEXT(PRED(P1)) = P2 
PRED(P2) = PRED(PI) 

P2 = NEXT(P2) 

END IF 

RETURN 

END 
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LOGICAL FUNCTION QO(P1, P2, START2, HEAD, NEXT) 

C ******ttt*****tt***t*******t***t*********************t****t********t** 

C LOGICAL FUNCTION TO DETERMINE WHEN THE END OF EITHER LIST IS REACHED 

C INPUT: POINTERS INTO TWO ALPHA SUBSTRINGS 

OUTPUT: LOGICAL VARIABLE WHICH IS TRUE IF EITHER POINTER IS AT THE 
END OF A STRING, AND FALSE OTHERWISE 

INTEGER HEAD(IOOO), NEXT(IOOO), PI, P2, START2 

GO = .TRUE. 

IF(P1.EQ.START2)THEN 
GO = .FALSE. 

RETURN 
END IF 

IF(NEXT(P2).EQ.0)THEN 
GO = .FALSE. 

RETURN 
END IF 

IF(HEAD(P2).GT.0)THEN 
GO = .FALSE. 

RETURN 
END IF 

RETURN 
END 
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SUBROUTINE MERGER (HEAD, NEXT, PI, P2, TAIL, PRED) 

C ♦♦♦*♦♦♦*♦*♦♦*♦♦♦♦*♦♦**♦♦♦**♦*♦♦*♦♦*♦*♦*♦♦*♦♦♦**♦♦♦♦♦♦♦♦♦*♦*♦♦***♦♦**♦♦ 

C THIS SUBROUTINE COMPARES TWO ELEMENTS - ONE IN EACH ALPHA SUBSTRING 
C AND PLACES THE ONE THAT IS THE SMALLEST ON THE LIST REPRESENTING THE 
0 NEW SORTED LIST 

C INPUT: HEAD, NEXT, PRED, AND POINTERS INTO THE h€AD ARRAY 

OUTPUT: HEAD, NEXT, PRED, AND POINTERS FOR THE NEXT TWO ELEMENTS 
TO BE CHECKED 

INTEGER HEAD(IOOO), NEXT(IOOO), PI, P2, TAIL , PRED(IOOO) 

IF(HEAD(P1).GT.HEAD(P2)) THEN 

CALL TACON1(NEXT, PI, TAIL, PRED) 

ELSE 

CALL TACON2(NEXT, P2, TAIL, PRED) 

END IF 

RETURN 
END 



SUBROUTINE TAOON1 (NEXT, PI, TAIL, PRED) 

Q !(c***:)c:)c***!)c*******!t::)c***:t:!(c****************************:)c:)c***********:(c*:t::)c* 

C ADDS THE ELEMENT TO WHICH PI POINTS TO THE TAIL OF THE NEW STRING 

C INPUT: NEXT, PRED, PI, TAIL 

C OUTPUT: NEXT, PRED, PI, TAIL 

INTEGER NEXT(IOOO), PI, TAIL, PRED(IOOO) 

NEXT(TAIL) = PI 
PRED(PI) = TAIL 
TAIL = PI 
PI = NEXT(PI) 

RETURN 

END 



90 



SUBROUTINE TAOON2 (NEXT, P2, TAIL, PRED) 

C ADDS THE ELEMENT TO WHICH P2 POINTS TO THE TAIL OF THE NEW STRING 

C INPUT: NEXT, PRED, P2, TAIL 

C OUTPUT: NEXT, PRED, P2, TAIL 

INTEGER NEXT(IOOO), P2, TAIL, PRED(IOOO) 

NEXT(TAIL) = P2 
PRED(P2) = TAIL 
TAIL = P2 
P2 = NEXT(P2) 

RETURN 

END 
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SUBROUTINE FXTAIL(HEAD, NEXT, PI, P2, TAIL, START2, PRED) 



C ATTACHES THE END OF THE LONGER SORTED STRING TO THE END OF THE MERGED 
C STRING 

C INPUT: PI, P2, HEAD, NEXT, PRED, TAIL 

C OUTPUT: PI, P2, HEAD, NEXT, PRED, TAIL, START2 

INTEGER HEAD(IOOO), NEXT(IOOO), PI, P2, TAIL, START2, P11, 

+ PRED (1000) 

IF(P1.EQ.START2)THEN 

CALL TAOON2(NEXT, P2, TAIL, PRED) 

RETURN 
END IF 

IF(NEXT(P2).EQ.0)THEN 

CALL TAOON1(NEXT, PI, TAIL, PRED) 

P11 = TAIL 

2 IF(P1 .NE.START2)THEN 

P11 = PI 
PI = NEXT(PI) 

GOTO 2 
END IF 

NEXT(PII) = P2 
PRED(P2) = P11 
RETURN 
END IF 

IF(HEAD(P2).GT.0)THEN 

CALL TAOON1(NEXT, PI, TAIL, PRED) 

P11 = TAIL 

3 IF(P1 .NE.START2) THEN 

P11 = PI 
PI = NEXT(PI) 

GOTO 3 
ENDIF 

NEXT(PII) = P2 
PRED(P2) = P11 
RETURN 
ENDIF 

RETURN 

END 
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SUBROUTINE CTREX(FL, NL, PL, FR, NR, PR, BSEED, PSEED) 

Q ♦*♦****♦****♦***♦♦♦**♦***♦****♦*********♦**♦♦♦♦♦♦**♦♦♦******♦♦**♦♦♦*♦♦ 

C THIS SUBROUTINE CHECKS F(C(R(P))) AGAINST K(F(R(P))) FOR EQUALITY. 

C INPUT: FUNDAMENTAL GROUP REPRESENTATIONS FOR BOTH ALGORITHMS 

OUTPUT: MESSAGE INDICATING THE OCCURRENCE OF UNEQUAL FUNDAMENTAL 

GROUPS 

INTEGER FL(IOOO), NL(IOOO), PL(IOOO), FR(IOOO), NR(IOOO), PR(IOOO) 

IF(NL(1).EQ.O.AND.NR(1).EQ.O) GO TO 11 

I = NL(1) 

J = NR(1) 

10 IF ( I .EQ.O.AND. J.EQ.O) GO TO 11 

IF (FL(O.EQ.O.AND.FR(J).EQ.O) GO TO 11 

IF (FL( I).EQ.FR(J)) THEN 
I = NL(I) 

J = NR(J) 

GOTO 10 

ELSE 

PRINT*, 'STRINGS UNEQUAL* 

PRINT*,’ ' 

IF(NL(1).EQ.0)THEN 

PRINT*, 'THE GROUP IS EMPTY AFTER ALG1 SEQUENCE’ 

ELSE 

PRINT*, ’THE FUNDAMENTAL GROUP AFTER ALG1 IS:’ 

CALL PRINTS(FL,NL) 

END IF 

IF(NR(1).EQ.O) THEN 

PRINT*, ’THE GROUP IS EMPTY AFTER ALG2 SEQUENCE’ 

ELSE 

PRINT*, ’THE FUNDAMENTAL GROUP AFTER ALG2 IS:’ 

CALL PRINTS(FR,NR) 

END IF 
PRINT*,’ ’ 

PRINT*, ’BSEED AND PSEED ARE’ , BSEED, PSEED 
END IF 

11 RETURN 
END 
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SUBROUTINE PR I NTS (HEAD, NEXT) 

C THIS SUBROUTINE CAN BE USED TO PRINT OUT ANY STRING THAT IS STORED IN 
C LINK LIST FORM 

C INPUT: LINK LIST ARRAYS 

C OUTPUT: HORIZONTAL STRING OF THE ELEMENTS ACCORDING TO THE NEXT() 

INTEGER HEAD(IOOO), NEXT(IOOO), STRING(IOOO) 

NSTART = NEXT(1) 

I = 0 

98 IF(NSTART.NE.O) THEN 
1 = 1 + 1 

STRING(I) = HEAD( NSTART) 

NSTART = NEXT (NSTART) 

GOTO 98 
END IF 

PRINT 111, (STRING(J), J = 1, I) 

111 FORMATC ■ ,2014) 

PRINT*,' ’ 

RETURN 

END 
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C ALL SUBRCUTINES BELOW THIS POINT WERE USED AS DEBUG TOOLS. THEY 
C GENERATE CRUDE GRAPHS WHICH PLOT THE BOARDS AND THE PATTIS. 

C TWE VALUES AQUIRED BY THE RAWSTRING , MERGE, AND MANY OTHER 
C SUBROUTINES WERE CHECKED USING THESE ROUTINES. 

SUBROUTINE GRAPH(BX, BY, X, Y, NOBS, NUMPTS) 

REAL*8 X(1000), Y(1000), BX(IOOO), BY(IOOO) 

CHARACTER*! MATR(53, 105) 

CALL GRAPPA (X, Y, NUMPTS, MATR) 

CALL GRAFOB(BX, BY, NOBS, MATR) 

DO 1 I = 1,53 

1 PRINT 111, (MATR(I,J), J=1,105) 

111 FORMATC ’,105A1) 

RETURN 

END 
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SUBROUTINE GRAPPA (X, Y, NUMPTS, MATR) 

Q itit******************************************^****^***^^**************** 

REAL+8 X(1000), Y(1000) 

CHARACTER*! MATR (53, 105), CH 

DO 10 L = 2,52 

DO 10 M = 2,104 
10 MATR(L,M) = ’ • 

DO 1 1=1, 105 

MATRd , I ) = ’*' 

1 MATR(53, I ) = **' 

DO 4 I = 1 , 53 

MATR(I ,1) = '*' 

4 MATR( I , 105) = '*' 

XK = X(1) 

YK = Y(1) 

CALL OOOR(XK, YK, IX, lY) 

CALL CHARPA(1, CH) 

MATR(IY,IX) = CH 

DO 2 K = 2, NUMPTS 

2 CALL FILL(K, X, Y, MATR) 

DO 3 K = 2, NUMPTS 
XK = X(K) 

YK = Y(K) 

CALL COOR(XK, YK, IX, lY) 

CALL CHARPA(K, CH) 

3 MATR(IY, IX) = CH 

RETURN 

END 



SUBROUTINE OOOR(BX, BY, IX, lY) 

REAL*8 BX, BY 

IX = MESHP(BX) * 2 
lY = 54 - MESHP(BY) 

RETURN 

END 
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INTEGER FUNCTION MESHP(X) 
REAL*8 X 

MESHP = ((X+1.)*(51./2.)) + 2. 

RETURN 

END 



SUBROUTINE CHARPA(K, CH) 

CHARACTER*! CH 

IF(K.EQ.I) THEN 
CH = *A’ 

RETURN 

END IF 

IF(K.EQ.2) THEN 
CH = ’B’ 

RETURN 

END IF 

IF(K.EQ.3) THEN 
CH = ’C’ 

RETURN 

END IF 

IF(K.EQ.4) THEN 
CH = ’D’ 

RETURN 

END IF 

IF(K.EQ.5) THEN 
CH = 'E' 

RETURN 

END IF 

IF(K.EQ.6) THEN 
CH = ’A’ 

RETURN 

END IF 

END 
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SUBROUTINE FILL(K, X, Y, MATR) 

Q ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦*******************t****j*:t********:|c******************* 

REAL*8 X(1000), Y(1000), U, V, DX, DY, XK, YK, XKM1, YKM1 
CHARACTER*! MATR(53,105) 

XK = X(K) 

YK = Y(K) 

XKM1 = X(K-1) 

YKM1 = Y(K-1) 

CALL COOR(XK, YK, KX, KY) 

CALL COOR(XKM1, YKM1, KXM1 , KYM1) 

L = IABS(KX - KXM1) - 1 
M = IABS(KY - KYM1) - 1 
IF (M.GT.L) L = M 
DX = (XK - XKM1)/(L+1) 

DY = (YK - YKM1)/(L+1) 

DO 1 J = 1,L 

U = XKM1 + J*DX 
V = YKM1 + J*DY 
CALL COOR(U,V, lU, IV) 

1 MATR(IV, lU) = '+’ 

RETURN 

END 



SUBROUTINE GRAFOB (X, Y, NOBS, MATR) 

Q ♦♦♦******************************************************************* 

REAL*8 X(1000), Y(1000) 

CHARACTER*! MATR(53, !05), CH 

MATR(27, 53) = ’*’ 

DO 2 K = !,NOBS 
XK = X(K) 

YK = Y(K) 

CALL OOOR(XK, YK, IX, lY) 

CALL CHAR(K, CH) 

2 MATR(IY, IX) = CH 

RETURN 

END 
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SUBROUTINE CHAR(K, CH) 

CHARACTER*! CH 

IF(K.EQ.I) THEN 
CH = ' 1 ’ 

RETURN 
END IF 

IF(K.EQ.2) THEN 
CH = '2' 

RETURN 
END IF 

IF(K.EQ.3) THEN 
CH = ’3' 

RETURN 
END IF 

IF(K.EQ.4) THEN 
CH = '4' 

RETURN 
END IF 

IF(K.EQ.5) THEN 
CH = *5' 

RETURN 
END IF 

IF(K.EQ.6) THEN 
CH = '6' 

RETURN 
END IF 

IF(K.EQ.7) THEN 
CH = *7' 

RETURN 
END IF 

IF(K.EQ.8) THEN 
CH = ‘8’ 

RETURN 
END IF 

IF(K.EQ.9) THEN 
CH = '9' 

RETURN 
END IF 

IF(K.EQ.IO) THEN 
CH = 

RETURN 
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END IF 



IF(K.EQ.II) THEN 
CH = ’#■ 
RETURN 
END IF 

IF(K.EQ.12) THEN 
CH = ’$' 
RETURN 
END IF 

IF(K.EQ.13) THEN 
CH = 

RETURN 
END IF 

IF(K.EQ.14) THEN 
CH = '&■ 
RETURN 
END IF 

IF(K.EQ.IS) THEN 
CH = '<t' 
RETURN 
END IF 

IF(K.EQ.ie) THEN 
CH = V 
RETURN 
END IF 

IF(K.EQ. 17) THEN 
CH = A' 
RETURN 
END IF 

IF(K.EQ.18) THEN 
CH = ’<’ 
RETURN 
END IF 

IF(K.EQ.19) THEN 
CH = ’>’ 
RETURN 
END IF 

IF(K.EQ.20) THEN 
CH = ’?' 
RETURN 
END IF 
END 
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APPENDIX B. EXAMPLE: GENERATING CLASS NAMES 



This appendix illustrates the procedures used in 
traversing a ring graph in breadth first order to produce the 
names of candidate homotopy classes. The traversal is done 
following the rules set forth in Chapter III. 

The problem considered is that of determining the names of 
all candidate classes of paths whose shortest representative 
does not cross itself between points a and z in the plane 
containing two fixed obstacles. Figure B.l illustrates the 
reference frame for the region and the location of the points 
a and z . 




2 



«1 



Figure B.l Region with n = 2 
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Once the reference frame is established, the region is 
modeled with a ring graph as shown in Figure B.2. This graph 
is traversed beginning at node III to produce a list of all 
character strings representing class names and the vertex in 
which a path in that class ends. 



I 




Figure B.2 Ring Graph for n = 2 

For illustration purposes, as the algorithm is 
implemented, all k+1 long offspring of a given k long string 
will be listed. If a k+1 long string violates any of the 
"sieving" criteria presented in Chapter III, excluding 
consideration of the ending node, it will be marked with a (*) 
and the rule it violates will be stated. That string will 
then be withdrawn from further consideration in the algorithm. 
In addition to the criteria set forth in Chapter III, those 
strings which duplicate any previous string will also be 
labelled and discarded. 
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Upon termination of the algorithm, those strings not 
labelled with a (*) will be searched for strings which do not 
terminate at the desired node, node I. These will be labelled 
with an (n) . Those strings not labelled with a (*) or (n) 
upon completion of this final check will denote the desired 
list of class names. 

The list generated by the algorithm follows. 



Violation 

Code 


End 

Node 


String 


Violation 


n 


III 


e 




n 


II 






n 


IV 


^1 






I 


“ 2^1 






I 


02^1 




★ 


I 


0^02 


duplicates 1 




I 






n 


II 


0:2132 




n 


IV 


02132 132 




n 


IV 


1 — • 




n 


II 






★ 


IV 


O2O2 O2 


cancel latio: 


n 


IV 


O2O2 ^2 




n 


IV 


02132 Cfc 2 




n 


II 






n 


II 


02^2 




★ 


II 


02132 


wraps 


n 


III 


02^2 ^ 1^2 




★ 


I 


02^2 ^ 2^2 


wraps 


n 


III 


O2B2 ^ 2^1 






I 


02132 02^2 
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★ 


III 


Q:2B^ ^2^1 


dupl i cates a26202a 


★ 


I 


ct2^2 


wraps 


★ 


III 


0:20:2 fi|02 


dupl icates a262O20: 


n 


III 


0:20:^ 62^1 




★ 


I 


0:2Oj fi2^2 


wraps 


★ 


III 


0:262 ^2^1 


dupl icates a20^iSi2^ 


★ 


I 


0:2^2 ^2^2 


wraps 




I 


0262^1®! 




★ 


III 


0:2^2 ^2^2 


dupl icates a2a2620 


★ 


I 


0:2^2 


wraps 


★ 


III 


0:2^2 ®1^2 


duplicates a26262a 


★ 


IV 


0262 ^2^2 ^1 


cancel lation 


★ 


II 


0:2^2 ^2^1 ^2 


wraps 


★ 


IV 


0:2^2 ^2^2 ^2 


wraps 


* 


1 1 


02^2 0:2^2 0:2 


crosses 


n 


II 


0:2^2 ^2^2 ^1 




* 


1 1 


0:20: 2 ^2^1 ^2 


wraps 


★ 


II 


0:262 0:262 0:2 


wraps 


* 


IV 


0:262 ^2^1 ^2 


crosses 


n 


IV 


0:262 0:262 62 




★ 


I 


0 2 62 0:262 ®2^1 


crosses 


★ 


III 


0:262 0:262 ^2^2 


wraps 


★ 


III 


0:262 0:262 ®2^1 


wraps 


★ 


I 


0:262 ^2^1 ^2^2 


crosses 



From the above list, the remaining candidate homotopy 
class names are: 



1. 




2. 


Qj^l 


3. 


Qjflj 


4. 


0: 262 0: 262 


5. 


0:262 0: 262 
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Figures B.3 through B.7 illustrate paths which belong to 



each of these classes. 




2 



a 



1 



Figure B.3 Class ct 2 B^ 
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z 



I 



p 



1 



p 




2 



a 



1 



Figure B.4 Class OjOj 



P 
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1 



Figure B.7 Class QjI32aj^fi^ 
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APPENDIX C. EXAMPLE; COUNTING NON- ISOMORPHIC GRAPHS 



The material of this appendix is presented to illustrate 
the ideas introduced in Chapter V. The example used counts 
the number of non-isomorphic ring graphs given n=7 using 
Burnside's Lemma and then demonstrates the procedures used in 
determining the number 'of ring graphs which may be generated 
given a particular partition of n. 

Given the number of obstacles n = 7, a regular 14-gon, as 
shown in Figure C.l is produced. Edges are numbered clockwise 
for convenience. 




Figure C.l 14-gon resulting from n=7 

The edges of the polygon are then colored black (bold) or 
white (fine) subject to the restriction that opposite edges 
are colored differently. The question to be answered is this: 
How many non-equivalent two-colorings exist? 
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The 4n 



28 total permissible permutations consist of 



2n = 14 rotations of the polygon, n = 7 reflections about 

diagonals drawn through opposite vertices, and n = 7 

reflections about bisectors drawn through opposite edges. The 
cycle notation for each of the 28 permutations is shown below 
along with the number of colorings left fixed by that 
permutation. These numbers are obtained by applying the 
formulas introduced in Chapter V. In cases where no coloring 
is left fixed by a permutation with more than one cycle, 
cycles containing opposite edges will be underlined. 



Rotations 

X, : (1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12){13)(14) - f(xo) = 2’ 

Xi : (1 2 3 4 5 6 7 8 9 10 11 12 13 14) • T{xi) = 0 

X, : (1 3 5 7 9 11 13)(2 4 6 8 10 12 14) - T(xj) = 2 

X, : (1 4 7 10 13 2 5 8 11 14 3 6 9 12) • f(x,) = 0 

X, : (1 5 9 13 3 7 11)(2 6 10 14 4 8 12) - T{x<) = 2 

Xj : (1 6 11 2 7 12 3 8 13 4 9 14 5 10) - T(Xs) = 0 

x« : (1 7 13 5 11 3 9){2 8 14 6 12 4 10) ■* T{Xj) = 2 

X, : (1 8)(2 9)(3 10)(4 11)(5 12)(6 13)(7 14) - T(x,) = 0 

X, : (1 9 3 11 5 13 7)(2 10 4 12 6 14 8) - f(x,) = 2 

X, : (1 10 5 14 9 4 13 8 3 12 7 2 11 6) - f{x,) = 0 

Xi, : (1 11 7 3 13 9 5)(2 12 8 4 14 10 6) • T(xi,) = 2 

Xii : (1 12 9 6 3 14 11 8 5 2 13 10 7 4) ■* T(Xu) = 0 

Xij : (1 13 11 9 7 5 3)(2 14 12 10 8 6 4) - f(xu) = 2 

Xij : (1 14 13 12 11 10 9 8 7 6 5 4 3 2) - T(xu) = 0 
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Reflec±ions About Diagcaials 

TCn : (1 14)(2 13)(3 12) (4 11) (5 10)(6 9)(7 8) - = 0 
7115 : a 2)(3 14)(4 13) (5 12) (6 11)(7 10)(8 9) • = 0 
7 Ci 5 : (2 3)(1 4)(5 14) (6 13) (7 12)(8 11)(9 10) • T(7Cu) = 0 
TCi, : (3 4)(2 5)(1 6) (7 14) (8 13)(9 12)(10 11) - f(7ii,) = 0 
Tti, : (4 5)(3 6)(2 7) (1 8) (9 14)(10 13)(11 12) • f(7Ci,) = 0 
Til, ; (5 6)(4 7)(3 8) (2 9) (1 10)(11 14)(12 13) • T(tCi,) = 0 

71,, : (6 7)(5 8)(4 9) (3 10) (2 11)(1 12)(13 14) - f(7tj,) = 0 



Reflections About Bisectors 

71,1 : (1)(8)(2 14)(3 13)(4 12)(5 11)(6 10)(7 9) • f(7i,i) = 2‘ 

71,, : (2)(9)(1 3)(4 14)(5 13)(6 12)(7 11)(8 10) - T(ti„) = 2* 

71,3 : (3)(10)(2 4)(1 5)(6 14){7 13)(8 12)(9 11) - fCTt,,) = 2< 

71,5 : (4)(11)(3 5)(2 6)(1 7)(8 14)(9 13)(10 12) • f(7i„) = 2* 

7C,5 : (5)(12)(4 6)(3 7)(2 8)(1 9)(10 14)(11 13) - f(7t,5) = 2< 

7t,5 : (6)(13)(5 7)(4 8){3 9)(2 10)(1 11)(12 14) - f(7C,5) = 2' 

71,, : (7)(14)(6 8)(5 9)(4 10)(3 11)(2 12)(1 13) • f(7t„) = 2‘ 

Applying Burnside's Lemma to the above permutations gives 
N = (1/28) [2U2 + 2 + 2 + 2 + 2 + 2 + 2^ + 2^+2^ + 2^2^ + 2^+2^] 
or, N = 9. 

Thus, the number of non-equivalent edge two-colorings of 
the polygon, and hence the number of non-isomorphic ring 
graphs for n = 7 fixed obstacles is obtained. 
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Next, we consider the partitions of the integer 7. They 
are given below with an arrow (-) indicating those with an odd 
number of summands. 



- 7 
6 1 
5 2 

-511 
4 3 

-421 
4 111 
-331 
-322 
3 2 11 
-31111 
2 2 2 1 
-22111 
2 11111 
-1111111 



Consider the partition 2, 2, 1,1,1 with k = 5 summands. 
Applying the information from Chapter V, a pentagon (since 
k = 5) is constructed as shown in Figure C.2. Vertices are 
numbered clockwise for convenience. 




Figure C.2 Pentagon Corresponding to Partition where k = 5 
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The 2k = 10 permutations of the pentagon consist of five 
rotations and five reflections about axes drawn through a 
vertex and bisecting its opposite edge. Each permutation is 
listed below, along with its cycle notation and cycle 
structure representation. 



Rotations 

X, : (1)(2)(3)(4)(5) - jq’ 

7Ci : (1 2 3 4 5) • Xj* 

Xj : (1 3 5 2 4) • X 5 * 

Xj : (1 4 2 5 3) ■* Xj* 

X, ; (1 5 4 3 2) ■* X 5 * 



Reflections 




Xj : 


(1)(2 5)(3 4) 


- 


Xj : 


(1 3)(2)(4 5) 


- 


X, : 


(1 5)(2 4)(3) 


- 


X, : 


(1 2)(3 5)(4) 


- 


X, : 


(1 4)(2 3)(5) 


- 



These permutations result in the cycle index 

= is + 5x^x|) . 



In applying Polya's Theorem to determine the number of 
non-isomorphic ring graphs generated by this partition, the 
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number of colors used is two, since the partition consists 
only of the summands 2 and 1. Letting these colors be b and 
w for black and white respectively, 

P(.((b+w), (b^+w^),..., (b^+w^)) 

yields 

+ 4 + 5 (f)+v) (jb^+vt^) . 

From this index, the coefficient of b w is desired since the 
summand 2 appears twice in the partition and the summand 1 
appears three times. 

Expanding the index gives 

— [lOjb® + 10b*w + 20jb^w^ + 20b^v^ + lObw* + lOn^] . 

10 

Hence, the coefficient of bw is 2, which is. the number of 
non-equivalent two-colorings which may be produced from the 
partition 2, 2, 1,1,1. The two corresponding non-isomorphic 
ring graphs are shown in Figure C.3. This same procedure is 
applied to the remaining candidate partitions to find the 
number of non-isomorphic ring graphs which may be produced 
from each one. 
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Figure C.3 Non-isoiiiorphic Ring Graphs Produced Using 
Summands 2, 2, 1,1,1 
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